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ABSTRACT 
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propagation,  the  method  uses  five  receiving  antennas  in  a  ladder  configuration.  Ratios 

of  received  signal  powers  are  compared  with  values  from  lookup  tables  to  determine 

the  correct  location  of  the  transmitter.  Computer  simulation  results  are  presented,  to 
verify  the  suggested  method. 
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I.  INTRODUCTION 


The  improved  characteristics  of  a  modern  missile  make  the  interception  process 
difficult:  the  reduced  Radar  Cross  Section  (RCS)  makes  it  hard  to  detect  by  using 
a  radar,  while  late  emission  leaves  little  time  for  intercepting  and  initiating  any 
necessary  defense  and  Electronic  Counter  Measures  (ECM)  activities. 

Conventional  interception  methods  give  only  direction,  and  rough  estimated 
range,  and  require  a  relatively  long  period  of  time. 

This  work  suggests  a  unique,  non-conventional  method  for  passive  range  and 
height  estimation  of  a  cruising  missile,  as  well  as  any  other  microwave  transmitter. 
The  method  is  based  on  over  sea  multipath  propagation,  and  uses  five  receiving 
antennas  in  a  ladder  configuration.  For  each  receiving  antenna  there  is  a  particu¬ 
lar  multipath  behavior,  different  than  the  other  antennas.  Comparing  theoretical 
ratios  of  antenna  pairs  with  actual  received  power  ratios  leads  to  an  estimated 
location  (range  and  height)  of  the  transmitter. 

The  motivation  for  the  need  of  such  a  method  is  given  in  Chapter  II.  Chapter 
III  describes  the  theory  behind  multipath  propagation  of  electromagnetic  waves, 
over  different  kinds  of  sea  surface.  It  also  explains  in  detail  the  model  for  rough 
sea  which  we  used  in  this  work. 

Chapter  IV  gives  the  basic  theory  of  how  the  multipath  propagation  can  be 
used  for  range  estimation,  and  the  conceptual  estimation  sequence.  The  detailed 
process,  the  assumptions  being  used  during  the  research,  and  simulation  results 
axe  given  in  Chapter  V. 
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Additional  parameters  that  were  not  taken  into  account  in  the  process,  as  well 
as  those  only  mentioned,  are  considered  in  Chapter  VI,  along  with  their  influence 
on  the  experiment.  The  results  are  analyzed. 

Chapter  VII  explains  the  computer  program,  and  shows  the  concept  behind 
the  software  process. 

Chapter  VIII  summarizes  and  concludes  the  whole  work,  and  gives  some  sug¬ 
gestions  for  future  steps. 

The  rough  sea  model  which  we  used  in  this  work  is  the  “Ament  Model.”  Other 
familiar  models  are  described  in  Appendix  A.  Listing  of  the  computer  program  used 
is  given  in  Appendix  B. 
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n.  NATURE  OF  THE  PROBLEM 


A.  GENERAL 

In  the  modem  battlefield,  where  th „  electromagnetic  complexity  and  density, 
as  well  as  a  missile’s  speed  are  very  high,  the  probability  of  intercept  is  low,  so 
detecting  and  intercepting  missions  became  hard  to  achieve.  The  policy  is  to  start 
the  transmission  as  close  to  the  target  as  possible,  and  to  transmit  for  a  minimum 
amount  of  time  in  order  to  be  less  “visible.” 

On  the  other  hand,  the  radar  cross  section  (RCS)  of  ships  and  especially 
missiles  is  being  reduced,  making  it  difficult  for  a  radar  to  detect  them. 

Therefore,  we  see  that  intercepting  a  missile,  especially  an  ultrasonic  one,  is  a 
tough  mission  and  might  be  achieved  only  during  the  last  seconds,  which  might  be 
too  late. 

By  using  a  conventional  interceptor,  even  if  it  is  a  fast  one  (monopulse  re¬ 
ceiver),  only  a  direction  can  be  found.  Using  the  ship’s  radar  for  range  finding  can 
be  worthless,  because  of  the  small  RCS  of  the  missile,  and  at  the  same  time  will 
help  the  other  side  to  detect  the  source. 

For  these  reasons,  we  prefer  to  find  a  missile  (as  well  as  other  ships’)  location 
(direction  and  range)  by  using  high  speed  passive  equipment. 

Two  common  methods  for  solving  this  problem  axe  being  used,  and  are  listed 
below: 

1.  Triangulation 

This  method  requires  two  independent  interceptors  on  two  different  plat¬ 
forms  sharing  the  information  of  the  target  direction.  The  location  of  the  target 


can  be  found  on  a  map,  as  shown  in  Figure  2.1.  Since  each  receiving  antenna  has 
its  beam  width,  the  position  estimate  is  actually  within  a  rhombus,  and  may  not  be 

as  accurate  as  desired.  The  direction  errors  of  the  interceptors  also  add  a  particular 
level  of  inaccuracy  to  the  estimation  process. 

Other  disadvantages  of  this  method  are  the  need  for  two  ships  to  be  in- 
vdved,  with  all  the  logistics  and  synchronization,  and  obviously  the  time  frame 
needed  for  this  process,  which  make  it  a  poor  choice  for  missile  defense. 


Interceptor  #  1 


Figure  2.1:  Triangulation  Geometry. 
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Interceptor  #  2 


Auto- triangulation 


In  order  to  estimate  the  location  by  using  a  single  interceptor  (single  ship), 
ve  can  make  the  triangulation  by  making  sequential  measurements,  using  the  fact 


that  the  interceptor  itself  is  moving  and  changing  its  place.  Obviously,  this  method 
takes  a  long  time  for  processing,  and  therefore  is  good  only  for  stable  targets  (shore 
radar,  for  example). 

Actually,  we  cannot  find  the  range  of  a  moving  target,  especially  a  cruising 
missile,  by  using  conventional  methods,  and  so  we  need  to  have  a  special  way  for 
estimating  the  transmitter  location. 

The  method  suggested  in  this  work  uses  the  multipath  propagation  of  a 
transmitted  electromagnetic  signal  in  order  to  estimate  the  transmitter  range  from 
the  receiver,  while  the  direction  is  still  being  found  using  any  common  method  such 
as  rotating  antenna,  Instantaneous  Direction  Finder,  etc. 

The  theory  of  a  general  multipath  propagation  over  smooth  and  rough  sea 
is  described  in  the  next  chapter. 


HI.  THEORY  OF  MULTIPATH  PROPAGATION 
OF  ELECTROMAGNETIC  WAVES 


In  a  usual  radar  case,  the  multipath  propagation  generally  causes  tracking 
difficulties,  especially  over  a  smooth  sea,  which  causes  specular  reflections. 

In  this  chapter  we  will  review  the  multipath  propagation  theory,  divided  into 
smooth  and  rough  terrains. 

A.  GENERAL 

A  transmitted  signed  from  a  radar  antenna  spreads  as  it  is  leaving  the  antenna 
surface.  A  direct  path  exists  between  the  transmitter  (radar)  antenna  and  the 
target  (receiver  antenna  in  our  case),  as  well  as  an  indirect  path  along  which  the 
signal  transmitted  from  the  transmitter  antenna  hits  the  sea  surface  and  is  reflected 
to  the  receiver  antenna. 

The  signals  arriving  from  these  two  paths  create  interference  in  the  receiver, 
and  the  relative  phase  difference  dictates  if  it  is  a  creative  or  destructive  inter¬ 
ference.  The  path  difference  depends  on  the  direct  range  R,  transmitter  antenna 
height  ht,  and  receiver  antenna  height  hr.  Other  factors  have  some  influence  as 
well,  and  will  be  discussed  later  in  the  next  sections.  Figure  3.1  shows  the  geometry 
of  the  multipath  propagation  over  a  smooth  sea,  assuming  a  flat  earth. 

Since  the  sea  state  has  significant  contribution  to  the  interference,  a  basic 
question  that  can  be  asked  is  how  to  distinguish  between  a  smooth  and  rough 
surface  (sea  in  our  case).  The  answer  for  this  question  is  not  a  straight-forward 
one,  and,  as  a  matter  of  fact,  there  is  neither  a  single  answer  nor  a  defined  border 
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IMAGE 

ANTENNA  f- - R - -| 

Figure  3.1:  Simplified  Multipath  Geometry  Over  Smooth  Terrain  and 
Flat  Earth  Assumption  [From  Ref.  1], 

between  the  two  cases.  A  method  for  having  a  quantitative  solution  for  this  problem 
was  suggested  by  Lord  Rayleigh,  and  is  described  herein  [after  Ref.  2]. 

Before  starting  the  detailed  analysis,  we  will  give  some  definitions  which  are 
relevant  to  our  discussion,  referring  to  Figure  3.2: 


Figure  3.2:  Geometry  of  Irregular  Surface  and  Path  Difference  [From 
Ref.  3]. 
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•  Wave  Height  Ah:  The  vertical  distance  between  the  wave’s  trough  and  its 
crest. 

•  Wavelength  /:  The  horizontal  distance  between  two  adjacent  crests  (measured 
in  the  direction  of  wave  propagation). 

•  Wave  Period  T:  Time  interval  between  two  adjacent  crests  passing  the  same 
fixed  point  relative  to  the  ground. 

Consider  an  irregular  surface  with  irregularities  of  height  Ah.  If  two  parallel 
rays  of  the  same  wavelength  A  hit  the  surface  with  a  grazing  angle  ip,  as  illustrated 
in  Figure  3.2,  the  path  difference  between  them  will  be  A r,  where: 


Ar  =  2Ahsintp.  (3-1) 


The  phase  difference  A<p  can  vary  between  0  and  7r,  where  a  totally  smooth 
surface  will  occur  when  A <f>  =  0.  Zero  phase  can  be  achieved  when  there  are  no 
irregularities,  i.e.,  Ah  =  0,  or  when  the  grazing  angle  is  zero.  From  the  above,  we 
can  assume  a  smooth  surface  when: 


Ah 

A 


0  or  ip  0. 


(3.3) 


When  the  phase  difference  between  the  direct  and  reflected  signal  components 
is  equal  to  7r,  there  will  be  a  destructive  interference,  and  no  received  energy  at  all 
if  there  is  a  perfect  reflection  (|p|  =  1). 


8 


XI 


In  order  to  choose  a  point  where  the  surface  becomes  rough,  we  can  arbitrarily 
choose  any  value  between  A 4>  =  0  and  A <t>  =  it.  The  most  popular  value  is  half  the 
way  (A <t>  =  7t/2),  as  was  chosen  by  Rayleigh,  and  hence,  the  “Rayleigh  criterion” 
for  a  smooth  surface  is: 


A  h  < 


(3.4) 


8  sin  rl> 

Other  values  have  been  chosen  as  well  [Ref.  2, 3, 4, 5, 6]. 

From  observing  Equation  (3.2),  we  can  learn  that  the  phase  behavior  is  a 
periodic  one,  hence  the  interference  will  have  maxima  and  minima  in  sequence. 

A  typical  received  signal  in  the  presence  of  multipath  propagation  over  a 
smooth  sea  is  shown  in  Figure  3.3.  The  amplitude  decay  is  due  to  the  spheri¬ 
cal  spreading  of  the  received  signal  (l/R1). 


A  :  !  £h  I'CC’  ~  SEA  S 
lev  ANT  i-ClCHT  16M 


Figure  3.3:  Multipath  Propagation. 
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B.  SMOOTH  SEA 


Looking  at  the  multipath  propagation  more  carefully,  we  should  notice  that 
there  are  additional  phase  and  amplitude  changes  in  the  reflected  wave  due  to  the 
reflection  itself,  as  well  as  amplitude  difference  between  the  reflected  wave  and  the 
direct  one,  since  the  propagation  path  lengths  are  not  the  same. 

The  amplitude  difference  due  to  the  difference  in  path  lengths  is  very  small, 
and  hence  can  be  neglected,  since  the  paths  differ  only  by  a  short  distance  compared 
to  the  total  path  length.  [Ref.  7] 

The  phase  and  the  amplitude  changes  due  to  the  reflection  itself  depend  on 
the  grazing  angle,  frequency,  and  polarization. 

Figures  3.4  and  3.5  illustrate  the  magnitude  and  phase  (respectively)  for  verti¬ 
cal  polarization,  while  Figure  3.6  shows  the  magnitude  of  the  reflection  coefficient 
for  horizontal  polarization. 

Note  the  small  changes  in  the  value  of  the  reflection  coefficient  for  horizontal 
polarization,  which  means  that  amplitude  changes  due  to  reflection  can  be  neglected 
for  small  grazing  angles  (assuming  reflection  coefficient  magnitude  =  1  ).  The  phase 
difference  in  horizontal  polarization  can  be  considered  as  180°,  since  for  grazing 
angles  of  0°  to  90°  it  changes  only  up  to  4°.  [Ref.  7] 

Assuming  a  flat  earth  and  a  smooth  terrain  (reflection  coefficient  =  —  1),  the 
ratio  between  the  received  power  in  the  receiving  antenna  in  the  presence  of  mul¬ 
tipath,  compared  to  a  free  space  (for  small  grazing  angles)  is  r]2  [Ref.  3,9],  where: 

rj2  =  4sm2 - — ^ —  0),  (3.5) 

using  the  same  notation  as  in  Figure  3.1.  The  shape  of  r /2  is  shown  in  Figure  3.3 
for  a  typical  case. 
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Figure  3.4:  Reflection  Coefficient  in  Vertical  Polarization 
[From  Ref.  7]. 

In  order  to  be  more  precise,  we  need  to  take  into  account  the  fact  that  the 
earth  is  not  flat.  This  fact  will  cause  some  changes  in  the  geometry,  and  hence  will 
influence  the  signal  amplitude. 

Actually,  the  reflection  coefficient  depends  on  three  main  parameters: 

1)  Plane  smooth  surface  reflection  coefficient  p. 

2)  Earth  curvature,  D, which  reduces  the  reflection  (in  some  cases,  depending  on 
the  geometry). 

3)  Surface  roughness  R„. 

The  actual  reflection  coefficient,  hence,  is  [Ref.  6]: 

p'  =  pDRn.  (3.6) 
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Figure  3.5:  Phase  Difference  in  Vertical  Polarization  [From  Ref.  7]. 

The  value  of  Rn  is  between  0  and  1,  depending  on  the  terrain,  and  it  will  be 
discussed  in  detail  in  the  next  section. 

The  earth  curvature  factor  D  (Divergence  factor)  has  also  values  between  0 
and  1,  where  in  the  flat  earth  case  it  will  be  unity. 

Practically,  in  the  scope  of  this  work,  we  can  use  the  flat  earth  assumption, 
since  the  geometry  of  the  problem  (transmitter  height  less  than  about  100  meter, 
receiver  antenna  height  at  about  the  same  and  both  heights  are  much  smaller  than 
the  distance  in  between)  leads  to  small  grazing  angle  and  allows  the  assumption  of 
parallel  direct  and  reflected  rays  [Ref.  7]. 

Even  though  the  transmitter  is  considered  a  point  source,  the  reflections  from 
the  sea  come  from  a  particular  zone,  rather  than  in  a  single  point.  We  shall  look 
for  all  the  points  (loci),  from  which  the  reflected  signal  will  reach  the  receiving 
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Figure  3.6:  Coefficient  in  Horieontal  Polarizati 


Figure  3.7:  First  Fresnel  Zone  [From  Ref.  2]. 

antenna  with  the  same  phase  shift,  relative  to  the  direct  radiation,  as  illustrated 
in  Figure  3.7.  The  loci  obey  Equation  (3.7). 


Ri  +  R.2  —  Rd  =  S, 


(3.7) 


or  after  arranging: 

Ri  +  R2  —  Rd  +  S  =  Constant.  (3-8) 

Equation  (3.8)  describes  ellipsoids,  so  that  the  reflection  zones  are  actually 
ellipsoids  which  form  the  Fresnel  zones  (Fresnel  ellipsoids).  Near  the  wavelength 
that  has  the  most  interest  to  us,  3  cm,  (or  frequencies  around  10  Ghz),  the  area 
of  the  first  Fresnel  zone  is  about  2000  square  meters  [Ref.  2,3,6].  It  is  common  to 
take  the  whole  first  Fresnel  zone  as  the  reflection  area. 

C.  ROUGH  SEA 

The  smooth  sea  case  is  generally  rare,  while  some  kind  of  roughness  is  more 
common,  producing  complications  to  the  problem  by  adding  more  parameters, 
some  of  them  arbitrary.  Due  to  these  facts,  the  rough  sea  case  is  difficult  to 
analyze  and  predict. 

A  more  accurate  solution  leads  to  more  complications  and  difficulties,  so  in 
practice  we  try  to  simplify  and  make  assumptions  in  order  to  get  a  solvable  case. 
Still  we  need  to  be  very  careful  so  that  the  result  will  approximate  reality  even 
after  using  these  assumptions. 

In  general,  we  can  say  that  the  first  problem  is  to  find  a  correct  compromise 
between  excess  rigor  which  may  yield  no  practical  result,  and  too  many  simplifica¬ 
tions  which  may  lead  to  incorrect  results  and  hence  wrong  conclusions  [Ref.  3]. 

In  practice,  the  common  general  assumptions  which  take  place  in  the  rough 
sea  case  Eire: 

•  Scattering  elements  are  much  larger  than  the  propagated  wavelength. 
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•  Only  the  fax  field  is  calculated, 

•  Diffraction  and  reflections  between  waves  are  neglected, 

•  A  particular  model  for  the  roughness  shape  is  assumed  (depending  on  the 
case),  and  average  wave  height  is  assumed,  and, 

•  Constant  sea  state  (roughness)  is  assumed  for  the  whole  area  of  interest. 

1.  Reflection  types 

At  this  point  we  will  distinguish  between  two  kinds  of  reflections: 

a.  Specular  reflection 

Similar  to  the  smooth  surface  (as  in  a  mirror  case),  specular  reflec¬ 
tion  can  be  characterized  as  obeying  the  classical  optics  rules,  as  well  as  having 
a  coherent  phase.  The  direct  source  for  this  kind  of  reflection  is  the  first  Fresnel 
zone  (Fresnel  ellipse)  radiated  points.  The  phases  of  the  different  points  (from  the 
receiver  point  of  view)  are  approximately  equal,  hence  the  reflection  is  coherent, 
as  mentioned  previously.  The  specular  reflection  fluctuations  are  relatively  small 
in  amplitude. 

b.  Diffuse  reflection  (scattering) 

Caused  due  to  irregularities  in  the  surface,  scattering  cannot  be  pre¬ 
dicted  from  the  phase  and  directivity  points  of  view.  This  scattering  will  occur 
over  larger  areas  than  the  specular  reflection  (firs*’  Fresnel  zone)  and  is  incoherent 
in  phase.  The  diffuse  scattering  has  large  amplitude  fluctuations.  We  can  also  say 
that  the  smooth  sea  analysis  is  deterministic,  while  the  rough  sea  analysis  (both 
specular  and  diffuse  reflections)  is  stochastic  in  nature  [Ref.  5]. 

In  general,  there  will  be  both  reflection  components,  since  a  surface  is 
neither  totally  smooth,  nor  extremely  rough.  The  two  components  are  illustrated 
in  Figure  3.8  [Ref.  2,7,10]. 
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The  specular  and  diffuse  reflection  coefficient  components  will  be  de¬ 
noted  as  pa  and  pj,  respectively.  The  effective  reflection  coefficient  will  be  a  mul¬ 
tiplication  of  the  smooth  terrain  coefficient  po,  and  the  relevant  specular  or  diffuse 
coefficient.  Hence,  for  the  specular  component,  pa,  the  reflection  coefficient  will  be: 

P  =  Po  ■  P$,  (3-9) 

and  obviously  will  be  smaller  than  that  of  a  smooth  sea  po  (since  each  coefficient 
is  less  than  unity). 

The  diffuse  reflection  will  have  a  coefficient  of 

p  =  po-pd,  (3.10) 


again,  smellier  than  po- 


ROUGH  TERRAIN 


Figure  3.8:  Illustration  of  Diffuse  and  Specular  Reflections 
[From  Ref.  10]. 
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The  diffuse  reflection  coefficient,  due  to  its  random  nature,  is  found 
to  have  a  Rayleigh  distribution  probability  function  for  the  amplitude,  as  given  in 
Equation  (3.11),  while  the  phase  has  a  uniform  distribution  between  0  and  27r. 
The  diffuse  reflection  coefficient  amplitude  has  the  form: 

pd  =  V2 -T  ■  exp(jffa),  (3.11) 

where: 


crQ  is  the  RMS  wave  slope  (unitless),  as  shown  in  Figure  3.2, 

T  =  0.77  •  [1  —  exp(— 47r<5)]  -  exp(— 4.73  ■  6),  and, 
crk  ■  sin  0 
S  =  A  ’ 


(3.12) 

(3.13) 


where  cu,  is  the  standard  deviation  of  the  wave  height,  related  to  the  peak  to  trough 
height.  [Ref.  5] 

Assuming  a  smooth  terrain,  there  will  be  a  center  point  of  reflection. 
The  diffused  reflection  comes  from  the  whole  area  surrounding  this  point,  the  glis¬ 
tening  surface  [Ref.  2,11],  which  depends  on  the  observation  angle  from  the  trans¬ 
mitter  to  the  receiver  ( Et  in  Figure  3.1),  and  on  the  RMS  wave  slope  (cra). 

A  typical  RMS  slope,  aa,  is  about  0.2  radians,  and  it  depends  on  the 
wind  velocity.  The  amplitude  of  the  reflection  coefficient  has  an  extreme  of  0.4, 
which  rarely  occurs,  while  a  typical  value  is  0.35.  [Ref.  11] 

We  should  notice  that  for  a  flat  earth  assumption,  the  grazing  angle 
0  is  the  same  as  the  observation  single  Et  from  the  transmitter  to  the  receiver 
(see  Figure  3.1).  The  idea  of  the  different  reflection  areas  due  to  observation  angle 
Et,  and  sea  state  (roughness)  is  demonstrated  in  Figure  3.9.  In  this  Figure,  aQ 
represents  the  rough  surface  RMS  slope. 
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The  specular  reflection  coefficient  has  been  analyzed  by  many,  and 
there  are  quite  a  few  models  for  this  component.  Each  model  uses  different  as¬ 
sumptions  for  the  computations,  and  therefore  there  are  (in  some  cases)  different 
results.  The  conclusion  from  this  is  that  the  results  are  highly  dependent  on  the 
roughness  type  which  is  assumed. 
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Figure  3.9:  Reflection  Areas  in  Different  Roughness  and  Observation 
Angie,  Et :  (a)  Specular  reflection,  os  ~  l,p  =  p0»  (h)  Slight 
roughness,  (<7*/A)  •  sin Et  cz  0.06,  (c)  Rough  Surface  [From 
Ref.  11]. 


In  order  to  be  sis  close  to  the  experimental  data  with  a  generic  formula¬ 
tion  as  possible,  we  can,  in  general,  define  the  reflection  coefficient,  p ,  as  [Ref.  12]: 


P  = 


cos  <f>  d<j>, 


(3.14) 
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where: 


<j>  is  the  phase  of  each  particular  element  of  the  reflected  wave  in  relation  to  a 
median  plane  surface,  and, 

p(<f>)  is  the  probability  distribution  of  all  the  phases  <f>  along  the  path. 

2.  Profile  Types 

Using  Equation  (3.14),  we  find  that  the  common  profile  models  for  rough 
terrain  are  [Ref.  2,12]: 

a.  Sawtooth  profile 

The  phase  <f>  varies  uniformly  between  maximum  values  of  ±0m,  and 
the  probability  distribution  therefore  is: 

P (<f>)  =  «t-.  (3-15) 

while  the  reflection  coefficient  is: 


P  = 


sin  <f>m 
4>m 


=  sinc0m. 


(3.16) 


The  idea  here  is  that  the  wave  height  is  uniform,  but  wave  slopes  are  arbitrary. 

b.  Sinewave  profile 

The  phase  <f>  behaves  like  a  sine  wave  with  maximum  amplitude  of 
±0m.  In  this  case,  the  probability  distribution  is: 


1 

2n<f>m  cos  v ' 


where  v  is  an  arbitrary  parameter,  and: 


(3.17) 


P  —  ‘fo(^m)) 

where  Jo  is  the  Bessel  function  of  zero  order. 
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p  =  exp  -(<72/2).  (3.19) 

The  latter  profile  is  similar,  in  a  sense,  to  the  Ament  model,  which  will  be  discussed 
later  in  this  chapter.  The  three  profiles  are  drawn  in  Figure  3.10  [Ref.  12]. 


o  TT  2.TT  3T  AT 


Figure  3.10:  Reflection  Coefficients  of  Common  Models  for  Irregular 
Terrain  [From  Ref.  12]. 
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A  further  discussion  on  other  models  is  given  in  Appendix  A. 

The  most  general  and  the  most  popular  model  is  the  one  by  W.  S. 
Ament.*  This  model  assumes  a  normal  (gaussian)  distribution,  and  it  is  popular 
since  it  does  not  include  any  discontinuities  as  other  distributions  do,  but  still 
meets  the  actual  (experimental)  data.  For  these  reasons  it  seems  to  be  the  most 
realistic  model  [Ref.  2]. 

The  Ament  model  is  discussed  and  analyzed  in  many  books  and  arti¬ 
cles,  and  from  now  on  we  will  concentrate  especially  on  it  [Ref.  2,  3,  5,  6,  7,  9,  11, 
13,  14,  15]. 

D.  AMENT  MODEL  FOR  IRREGULAR  TERRAIN 

The  Ament  model  assumes  that  the  wave  height  is  normally  (gaussian)  dis¬ 
tributed,  and  it  does  not  take  into  account  the  shadowing  effect,  where  parts  of  the 
surface  which  are  hidden  by  higher  parts  are  not  being  illuminated  by  the  incident 
wave.  The  shadowing  effect  t Hikes  place  at  very  small  grazing  angles,  while  for 
higher  angles  it  has  no  effect.  The  geometry  is  shown  in  Figure  3.11.  [Ref.  13] 
Figure  3.12  shows  the  influence  of  the  shadowing  effect  on  the  wave  height 
standard  deviation,  related  to  the  grazing  angle  [Ref.  13]. 

The  specular  component  of  the  reflection  coefficient  in  this  model  is: 

p,  =  exp -2  f - - - )  ,  (3.20) 

where. 

ah  is  the  RMS  deviation  of  the  waves  height, 
ip  is  the  grazing  angle  (refer  to  Figure  3.1),  and, 

*  W.  S.  Ament,  “Toward  a  Theory  of  Reflection  by  a  Rough  Surface,” 
Proc.  I.R.E.,  voi.  41,  p.  142,  January  1953. 
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Figure  3.11:  Shadowing  Effect  Geometry  [From  Ref.  13], 


Figure  3.12:  Shadowing  Effect  on  Wave  Height  Standard  Deviation 
[From  Ref.  13]. 

A  is  the  incident  wavelength. 

Figure  3.13  shows  the  specular  reflection  coefficient  component  according  to 
Ament  model  [Ref.  7]. 

The  parameter  a is  proportional  to  the  sea  state,  and  is  given  with  other 
relevant  parameters  in  Table  3.1  [after  Ref.  2,11]. 


▼ 


Figure  3.13:  Rough  Surface  Reflection  Coefficient  According  to  Ament 
Model. 


TABLE  3.1:  WAVE  HEIGHT  PARAMETERS  AND  CRITICAL 
ANGLES  FOR  SPECULAR  REFLECTION,  VERSUS  SEA  STATE. 


sea  state 

number 

description 
of  sea 

wave 

height  (m) 

RMS  height 
Oh  (m) 

critical  anj 
3  GHz 

gle  (  ),  V’max 

10  GHz 

0 

calm 

0 

0 

1 

smooth 

0-0.3 

0-0.065 

>  6 

>  1.8 

2 

slight 

0.3-1 

0.065-0.21 

1.8-6 

0.5-1. 8 

3 

moderate 

1-1.5 

0.21-0.32 

1.2-1. 8 

0.3-0. 5 

4 

rough 

1. 5-2.5 

0.32-0.54 

0.7-1. 2 

0.2-0. 3 

5 

very  rough 

2.5-4 

0.54-0.86 

0.4-0. 7 

0.12-0.2 

6 

high 

4-6 

0.86-1.3 

0.3-0. 4 

0.09-0.12 

7 

[  very  high 

[  >  6 

|  >1.3 

<  0.3 

<  0.09 
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Using  the  Ament  model,  we  can  distinguish  more  carefully  between  rough  and 
smooth  sea,  and  check  the  conditions  for  diffuse  and  specular  reflection.  Recall 
(3.20)  for  convenience: 

„  /2tt crh  sin 

Pi  —  exp~2  ^  -  J  . 

We  can  consider  the  terrain  as  a  rough  one  if  p,  <  0.5,  which  is  equivalent  to: 

^  sin tAm ax  >  0.0665.  (3.21) 

For  a  particular  wavelength  at  a  given  sea  state,  there  is  a  critical  grazing 
angle  0max,  calculated  from  (3.21),  above  for  which  the  sea  surface  is  taken  as 
rough.  This  angle  is  shown  in  Table  3.1. 

The  specular  reflection  coefficient  component  is  significant  until  the  grazing 
angle  exceeds  twice  the  critical  angle.  Beyond  that  point,  the  diffuse  component 
becomes  the  significant  one.  The  different  regions  of  roughness  type  are  illustrated 
in  Figure  3.14.  [Ref.  2,11,14] 
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Another  way  of  looking  at  the  specular  (coherent)  and  diffuse  (incoherent) 
scattering  regions  is  given  in  Figure  3.15.  Here  we  can  see  the  boundaries  between 
the  scattering  regions  for  different  grazing  angles.  The  graph  in  Figure  3.15  is  given 
for  2  GHz.  For  higher  frequencies,  the  graph  will  become  more  compressed,  by  the 
reciprocal  of  the  frequency.  [Ref.  15] 


Figure  3.15:  Scattering  Regions  Envelopes  at  2  GHz  [From  Ref.  15]. 

A  similar  graph  for  a  constant  wave  height  of  0.92  meter  and  different  frequen¬ 
cies  is  given  in  Figure  3.16.  Here  again,  for  higher  wave  heights,  the  envelopes  will 
be  compressed  by  the  same  factor  by  which  the  height  is  increased.  [Ref.  15] 
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Figure  3.16:  Scattering  Regions  Envelopes  at  Peak  to  Through  Wave 
Height  of  0.92  meter  [From  Ref.  15]. 
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IV.  THEORY  OF  RANGE  ESTIMATION 
USING  MULTIPATH  EFFECT 


A.  THEORETICAL  SOLUTION 

The  multipath  effect  attenuation,  as  a  function  of  range  over  a  smooth  sea,  is 
well  known,  if  both  the  transmitter  and  receiver  height,  as  well  as  the  frequency  of 
the  signal  are  given.  This  was  described  in  the  previous  chapter. 

For  a  rough  sea,  the  overall  shape  remains  similar,  except  that  the  curve  has 
more  fluctuations.  Still,  by  using  a  particular  model,  such  as  the  Ament  model  in 
our  case,  we  can  deal  with  the  multipath  propagation  as  a  known  factor. 

Figures  4. 1  and  4.2  show  typical  multipath  propagation  over  a  calm  sea  (state 
0)  and  over  a  very  high  sea  (state  7),  respectively.  We  can  notice  that  in  the  rough 
sea  case,  the  attenuation  coefficient  decreases  arbitrarily  in  the  short  ranges,  and 
has  an  arbitrary  ripple  over  the  whole  range. 

From  the  above,  the  multipath  attenuation  rj 2  can  be  written  as  a  function  of 
the  range  between  the  transmitter  and  the  receiver  R  and  of  the  transmitter  height 
ht,  for  a  given  receiver  height  hr: 

rj2  =f(R,ht)\  .  (4.1) 

The  received  signal  is  attenuated  relative  to  free  space.  It  is  a  combination  of 
both  the  direct  and  indirect  components.  Hence,  the  received  signal  can  be  written 
as: 

Sr  =  Sd-TJ2,  (4.2) 
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where: 


Sr  is  the  received  signal  power, 

Sd  is  the  direct  component  of  the  signal  power,  and, 

r)2  is  the  multipath  power  attenuation  factor,  given  by  Equation  (3.5). 

Since,  as  mentioned  before,  we  cannot  specify  the  direct  component  Sd,  we 
need  to  eliminate  it  in  order  to  solve  the  problem.  A  simple  way  for  doing  that  is 
by  using  a  pair  of  receiving  antennas,  one  above  the  other,  so  that  the  range  R  will 
remain  the  same,  while  their  heights  have  no  practical  influence  on  the  direct  signal 
amplitude.  We  should  also  remember  that  the  transmitter  height,  ht,  is  constant 
for  both  receiving  antennas. 

If  we  use  the  ratio  of  the  output  powers  from  the  two  antennas, 

(Sr)i  _  sd-ni  _  nl 

(Srh  Sd  •  r\\  vV  1  ’ 

we  will  be  dependent  only  on  the  geometry  of  the  receiver  and  the  transmitter, 
without  caring  about  the  direct  signals’  amplitude. 

Referring  to  Equation  (4.1),  we  notice  that  we  still  have  two  unknowns,  R  and 
ht.  Using  the  fact  that  we  have  two  receiving  antennas,  while  the  range  R  and  the 
transmitter  height  ht  are  the  same,  we  need  to  find  the  correct  ratio  which 

will  give  the  right  combination  of  range  and  height  (R  and  ht).  Figure  4.3  shows 
the  geometry  of  the  problem.  A  more  detailed  explanation  of  the  sequence  is  given 
below. 

Let  us  assume,  for  simplicity,  that  the  multipath  propagation  curve  is  a  linear 
function  (we  shall  return  to  the  actual  case  later  on),  hence  it  behaves  like  an 
isomorphic  function. 
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Figure  4.1:  Multipath  Propagation  Over  a  Calm  Sea. 
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Figure  4.3:  Range  Estimation  Basic  Geometry. 

The  power  ratio  of  the  outputs  of  the  two  receiving  antennas  can  be  written 
as: 

(Sr)i_i7?  _Jhn\ 

(Srh  nl  f(R,ht)\hr  2  '{h*)' 

so  it  actually  depends  on  the  receiving  antenna  heights,  which  axe  known. 

Since  we  assumed  that  the  problem  has  an  isomorphic  characteristic,  a  given 

ratio  has  only  a  single  combination  of  range  and  transmitter  height  which  fits  this 

ratio. 

Keeping  in  mind  those  conclusions,  the  range  estimation  process  will  require: 

1 )  receiving  a  signal  and  analyzing  its  frequency  and  direction. 

2)  using  Equation  (3.5),  creating  two  (off  line)  tables  (matrices)  at  that  particular 
frequency,  for  which  the  coordinates  are  range  and  transmitter  height,  while 
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the  contents  axe  the  attenuation  coefficients  for  each  R  and  ht  combination, 
while  the  parameter  is  the  receiver  height  (hence  there  are  two  tables). 

3)  calculating  an  (off  line)  look-up  table  for  the  same  frequency,  ranges  and  trans¬ 
mitter  heights,  which  has  the  values  of  the  ratio  between  the  two  tables  created 
in  step  2. 

4)  calculating  the  ratio  between  the  two  received  signals,  which  is  similiar  (see 
Equation  4.3)  to  the  ratio  of  the  attenuation  coefficients. 

Now,  all  we  have  to  do  is  to  search  in  the  look  up  table  from  step  3,  and  find 
the  range  and  the  transmitter  height,  so  the  last  step  will  be: 

5)  searching  for  the  right  combination  of  range  and  height. 

B.  PRACTICAL  SOLUTION  (OVER  SMOOTH  SEA) 

Practically,  the  multipath  propagation  leads  to  a  multi-valued  curve,  hence 
the  problem  is  not  isomorphic  any  more,  but  has  an  ambiguous  solution.  Since 
the  function  is  not  linear  any  more,  and  involved  with  two  unknown  variables, 
we  cannot  use  any  analytical  solution  (like  we  could  have  used  in  the  theoretical 
case),  but  only  a  numerical  solution  can  take  place.  The  idea  is  similar  to  the  one 
described  in  the  theoretical  case  with  a  pair  of  antennas. 

We  saw  before  that  two  receiving  antennas  give  a  single  ratio,  hence  we  cannot 
solve  the  ambiguity  problem  (because  of  the  periodic  behavior  of  the  multipath 
effect). 

Adding  a  third  antenna  will  give  a  second  relation,  so  we  can  now  use  the 
process  described  previously  for  the  two  pairs,  with  a  few  candidate  combinations 
( R  and  ht),  and  by  comparing  their  locations  between  the  two  lists  of  candidates, 
choosing  that  one  which  appears  in  both  of  the  lists;  all  the  other  locations  are 
false  ones. 
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If  the  solution  still  is  not  a  single  one,  we  can  add  another  antenna(s),  until 
the  comparison  between  all  of  the  ratios  leads  to  a  single  location.  An  important 
condition  is  that  all  of  the  receiving  antennas  be  installed  one  above  the  other. 

A  trade-off  between  the  number  of  antennas  and  accuracy  should  be  done: 
on  one  hand  we  want  minimum  ambiguous  points  (different  points  with  the  same 
ratios  group)  which  dictates  many  antennas,  but  on  the  other  hand  we  should  use 
a  realistic  number  of  antennas. 

The  influence  of  the  different  parameters  on  the  multipath  behavior  can  be 
seen  in  Figures  4.4  to  4.6. 

From  observing  these  figures,  we  can  conclude  that: 

1)  as  the  frequency  increases,  the  curves  become  more  complex  (signal  strength 
varies  more  rapidly  with  range), 

2)  complexity  increases  with  transmitter  height,  and, 

3)  for  increased  receiving  antenna  height,  the  curves  become  more  complex.  This 
fact  suggests  keeping  the  antennas  at  lower  height,  but  on  the  other  hand,  we 
should  separate  them  enough  to  reduce  the  ambiguity.  The  receiving  antennas 
in  Figure  4.6  are  0.2  meters  apart,  so  the  curves  are  very  similar  to  each  other. 
In  Figure  4.7,  the  antennas  are  separated  by  3  meters,  and  the  curves  are  well 
distinguished. 


Note  that  for  short  range  we  have  very  high  complexity,  which  will  be  treated  later 


on. 


C.  PRACTICAL  SOLUTION  (OVER  ROUGH  SEA) 

In  the  case  of  rough  sea,  the  main  idea  remains  the  same,  although  the  signal 
exhibits  fading  as  a  function  of  the  sea  state.  As  the  receiving  antenna  is  raised, 
the  influence  of  the  roughness  becomes  more  important.  Figure  4.8  shows  the 
multipath  propagation  curves  over  sea  state  7,  for  different  antenna  heights.  It 


33 


is  obvious  from  these  graphs  that  the  higher  the  antenna  is,  the  more  rapid  the 
fading. 

The  look-up  tables  should  take  in  account  also  the  sea  roughness  parameter, 
where  in  our  case  we  choose  the  Ament  model,  as  described  in  Chapter  III. 

After  considering  the  facts  mentioned  above,  the  estimation  process  remains 
the  same  as  described  in  the  previous  sections. 
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Figure  4.4:  Variation  of  Signal  Strength  with  Range,  for  Different 
Transmitter  Heights. 
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Figure  4.5:  Variation  of  Signal  Strength  with  Range,  for  Different 
Frequencies. 
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Figure  4.7:  Variation  of  Sig 
Receiver  Height 
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Figure  4.8:  Variation  of  Signal  Strength  with  Range,  for  Very  High  Sea 
and  Different  Receiver  Antenna  Heights. 
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Now  that  we  have  presented  the  theory  of  how  to  use  the  multipath  effect 
for  estimating  range,  we  can  discuss  how  we  applied  that  for  solving  the  problem 
in  our  research.  The  assumptions,  practical  difficulties,  and  then  some  simulation 
results  will  be  presented. 


A.  ASSUMPTIONS 

Since  the  actual  problem  is  a  very  complicated  one,  due  to  many  surface 
parameters,  where  some  of  them  are  unknowns  and/or  very  difficult  to  predict,  we 
made  some  logical  assumptions  in  order  to  obtain  a  solution. 

•  Polarization:  the  seeker  was  assumed  to  have  horizontal  polarization,  hence 
the  reflection  coefficient  was  taken  as  unity  in  amplitude,  with  a  180°  phase 
shift.  The  differences  for  vertical  polarization  will  be  discussed  in  Chapter  VI. 

•  Sea  state:  the  whole  area  of  interest,  i.e.,  the  whole  range  between  the  trans¬ 
mitter  and  the  receiver  was  assumed  to  have  the  same  sea  state  condition. 

•  Earth  curvature:  a  plane  earth  was  assumed.  A  spherical  earth  case  is  dis¬ 
cussed  in  Chapter  VI. 

•  Atmosphere:  a  homogeneous  atmosphere  was  assumed. 

•  Movement:  the  suggested  method  considers  relative  movement  between  the 
transmitter  and  the  receiver.  For  homing  missile,  this  is  obvious.  In  other 
cases,  such  as  estimating  range  between  own  ship  and  a  stable  transmitter 
(shore  based  radar,  for  example),  we  shall  assume  that  the  ship  is  moving 
relative  to  the  radar.  If  there  is  no  relative  movement  at  all,  the  result  can  be 
incorrect. 

•  Geometry:  It  was  assumed  that  the  antenna  separation,  from  the  point  of 
reflection,  was  much  larger  than  the  signal  wavelength  [Ref.  2].  Since  we 
are  dealing  with  microwave  signals  (from  2  GHz  at  least),  this  assumption  is 
valid.  (If  this  condition  is  not  true,  the  reflection  coefficients  have  some  more 
complexity.) 
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•  Range  interval:  almost  any  range  interval  can  be  assumed,  but  note  that 
for  the  long  ranges  we  should  check  the  geometry  and  see  if  the  plane  earth 
assumption  is  still  applicable.  Since  we  concentrated  in  this  work  on  homing 
missiles,  the  seeker  transmitting  range  is  the  one  of  importance,  and  hence  a 
maximum  range  of  20  Km  (about  11  NM)  is  considered,  while  a  minimum  of 
1  Km  (about  half  a  NM)  assumed. 

B.  ESTIMATION  PROCESS 

With  the  assumptions  of  the  previous  section  in  mind,  let  us  repeat  briefly  the 
range  estimation  process.  We  shall  look  more  closely  at  the  details  shortly. 

1.  General  Discussion: 

A  general  solution  requires  the  following  steps: 

a.  Frequency  Measurement: 

Since  the  whole  process  is  frequency  dependent,  this  will  be  the  very 
first  step  made  by  any  EW  receiver. 

b.  Off  Line  Tables: 

In  this  step  we  build  the  theoretical  tables  for  all  ranges  and  heights  of 
interest  (i.e.,  the  possible  locations  of  the  transmitter),  for  each  receiving  antenna 
height.  The  inputs  for  creating  those  tables  are  frequency,  receiving  antenna  height, 
and  sea  state,  which  dictates  the  attenuation  factor  to  be  placed  in  the  table. 

c.  Attenuation  Coefficient  Ratio  Tables: 

After  we  have  the  theoretical  receiving  values  (from  the  previous  step), 
we  calculate  (for  each  range  and  height  combination)  the  ratios  of  the  output  signals 
for  each  receiving  antenna  pair. 

d.  Amplitude  Measurement: 

We  must  measure  the  received  signal  amplitude  of  each  receiving  an¬ 
tenna  (we  will  discuss  later  the  number  of  receiving  antennas),  and  make  a  table  of 
amplitude  ratios.  For  example:  suppose  we  have  three  receiving  antennas,  where 
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each  one  receives  a  signal  with  amplitude  Si,S2,  and  s3,  respectively.  The  ratio 
table  will  have  two  elements:  s 2/^1  and  s3/s\. 

e.  Comparison: 

Now,  that  we  have  the  measured  signal  ratio  (from  step  d),  and  the 
theoretical  one  (step  b),  we  need  to  compare  -  for  each  receiving  antenna  pair,  the 
actual  value  with  the  off  line  (theoretical)  value,  and  record  range-height  combina¬ 
tions  for  which  they  are  equal.  These  range-height  combinations  represent  possible 
transmitter  locations. 

Since,  as  we  mentioned  in  Chapter  IV,  the  problem  has  an  ambiguous 
solution,  there  will  be  several  possible  locations,  even  though  only  one  of  them  is 
the  correct  one.  The  next  step  is  to  determine  the  correct  location. 

f.  Validation: 

In  order  to  choose  the  correct  location  from  all  the  candidate  locations 
determined  in  step  e,  we  will  compare  all  the  locations  determined  for  the  first  pair 
of  receiving  antennas,  and  the  second  pair  of  receiving  antennas.  There  is,  ideally, 
only  a  single  location,  which  should  appear  common  to  all  the  pairs. 

Practically,  because  of  the  high  ambiguity  level  (especially  for  the  rel¬ 
atively  short  ranges),  and  also  because  of  approximations  in  the  implementation, 
there  will  still  be  more  than  a  single  location.  In  order  to  resolve  this  problem, 
we  shall  use  more  receiving  antennas,  which  provides  more  ratios.  This  allows  the 
ambiguity  to  be  reduced  to  an  acceptable  minimum. 

g.  Output: 

The  matrix  location  is  translated  into  an  estimate  of  transmitter  range 
and  height.  This  is  accomplished  by  using  the  coordinates  of  the  look  up  table. 
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2.  Detailed  Discussion 

The  following  section  explains  in  detail  every  step  listed  before. 

a.  Frequency  Measurement: 

Since  this  requires  an  external  system,  we  will  not  discuss  it.  Any  EW 
receiver  on  board  can  be  used  for  frequency  measurement  and  direction  finding. 
The  frequency  data  is  provided  as  an  input  to  the  range  estimation  process. 

b.  Off  Line  Tables: 

This  step  uses  the  theoretical  equations  for  calculating  the  multipath 
attenuation  factor  r]2 .  For  each  receiving  antenna  we  have  a  particular  table,  where 
the  coordinates  are  transmitter  height,  ht  (vertical  axis),  and  range,  R  (horizontal 
axis). 

The  height  axis  is  divided  into  36  “slices”  with  a  1  meter  resolution, 
starting  at  a  minimum  height  of  5  meters,  and  ending  at  a  maximum  of  40  meters. 
(These  numbers  are  parameters,  which  can  be  changed  easily.) 

The  range  axis  runs  from  1  Km  to  20.5  Km,  and  was  originally  divided 
into  130  equal  steps  of  50  meters  each.  Since  the  ambiguity  for  short  ranges  is 
very  high,  this  partition  caused  some  problems  in  finding  the  correct  location  for 
short  ranges.  This  suggested  the  need  for  higher  resolution  at  short  range.  The 
chosen  partition  is  130  steps  of  50  meters  each,  up  to  the  range  of  7.5  Km,  and 
then  reduced  resolution  of  125  steps,  100  meter  each,  until  the  maximum  range. 
(Again,  all  these  numbers  are  parameters  which  can  be  changed.) 

The  multipath  attenuation  factor  is  calculated  by  using  the  following 

formulae: 


The  basic  attenuation,  rj2 ,  for  a  smooth  sea  case  is  [Ref.  3,9]: 

n  .  .  9  2ltht(ht  “f”  ^r)  , 

V 2  =  4  sm2 - x!r—  (0-0), 


(5.1) 
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where:  ht  is  the  transmitting  antenna  height, 
hr  is  the  receiver  antenna  height, 

A  is  the  signal  wavelength, 

R  is  the  range  between  transmitter  and  receiver,  and, 
rp  is  the  grazing  angle. 

In  the  case  of  a  rough  sea,  we  have  to  multiply  77 2  by  the  reflection 
coefficient,  p ,  which  is  different  in  various  models.  In  our  case  we  used,  as  described 
in  Chapter  III,  the  Ament  model: 


p  =  exp  —2 


( 


sin  rp 
A 


(5.2) 


where:  <7/,  is  the  RMS  deviation  of  the  waves’  height, 
xp  is  the  grazing  angle,  and, 

A  is  the  signed  wavelength. 

Therefore,  in  a  rough  sea  case,  the  multipath  attenuation  factor  will 
be: 


trough  Vtmooth  '  P-  (5.3) 

sea  sea 


The  grazing  angle  rp  in  equation  (5.2)  is  a  derivative  of  the  geometry. 
It  depends  on  transmitter  height,  receiver  height  and  range.  Figure  5.1  shows  the 
geometry  from  which  we  calculate  rp. 

We  shall  now  derive  the  grazing  angle  rp  [According  to  Ref.  7].  From 
Figure  5.1  we  see  that 

Rd  =  y/(ht  ~  hr)2  +  R2.  (5.4) 


We  define  the  range  difference  between  the  indirect  path  and  the  direct 
path  as  6,  where: 


—  Ri  +  I?2  ~  Rd- 


(5.5) 
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Figure  5.1:  Flat  Earth  Multipath  Geometry  [From  Ref.  7]. 

Looking  at  the  image  of  the  transmit  antenna,  we  see  that 

Ri+R2  =  V(ht  +  hry  +  i?2.  (5.6) 


Substituting  Equations  (5.4)  and  (5.6)  into  Equation  (5.5),  we  obtain: 


6  =  y/( ht  +  hr)*  +  R 2  -  V(ht  -hr)2  +R2-  (5.7) 


Let  us  assume  that  the  range  R  is  much  greater  than  the  antenna 
heights,  and  also  R?  »  ( ht  +  hT)2 .  In  this  case,  we  can  approximate  Equations 
(5.4)  and  (5.6)  as: 


Rd  =  y/{ht  -  hr)2  +  R2  «  R 


1  _i_  +  ^r)2 

2R2 


(5.8) 
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Since  we  know  all  the  parameters  from  solving  previous  equations,  we 
can  now  find  rp.  (Actually,  we  do  not  know  the  range  R,  but  we  should  remember 
that  in  this  step  we  axe  building  tables,  so  the  range  is  being  used  as  a  parameter.) 

The  only  variable  still  missing  for  solving  Equation  (5.2),  (the  reflection 
coefficient),  is  the  wave  height  RMS  deviation  <7/,.  This  value  depends  on  the  sea 
state,  and  is  given  in  tables  in  the  literature.  Table  3.1  shows  this  parameter  (as 
well  as  some  others),  as  a  function  of  sea  state. 

The  propagation  attenuation  factor  (  in  the  one  way  case)  also 

has  its  impact  on  the  received  signal.  Since  we  deal  with  ratios,  this  factor  cancels 
itself,  and  therefore  had  not  been  taken  into  account  during  the  calculations. 
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c.  Attenuation  Coefficient  Ratio  Tables  (rjf/rjj): 

Having  created  the  different  tables  for  each  antenna,  we  now  create  a 
new  set  of  tables,  which  contain  the  ratios  of  the  received  signal  amplitudes  for 
the  antenna  pairs.  For  three  receiving  antennas  we  get  two  tables,  while  for  N 
receiving  antennas  there  will  be  N  —  1  tables. 

Figure  5.2  shows  graphs  of  the  ratio  table  data  for  sea  state  zero, 
at  9  GHz  and  a  transmitter  height  of  20  meters.  Note  that  for  clarity  only  one 
“slice”  is  shown,  referred  to  the  particular  transmitter  height.  The  table  is  actually 
three-dimensional  (range,  transmitter  height,  and  amplitude  ratio). 

d.  Amplitude  Measurement: 

The  first  and  major  assumption  is  that  all  the  receiving  channels  are 
calibrated.  In  this  case,  two  similar  signal  amplitudes  at  the  input  to  different 
channels  will  create  the  same  amplitude  at  the  outputs  of  those  channels. 

Actually,  there  is  no  need  for  single  channel  amplitude  measurements, 
since  we  deal  with  ratios.  For  this  reason,  we  measure  the  individual  signal  ampli¬ 
tudes  and  immediately  divide  between  pairs  to  form  the  desired  ratios. 

e.  Comparison: 

We  now  have  in  hand  all  the  details  needed  to  find  the  correct  location 
of  the  transmitter,  assuming  of  course,  that  the  model  is  accurate. 

A  comparison  is  made  among  the  ratio  levels  from  steps  c  and  d.  Every 
place  that  both  levels  are  the  same,  is  “marked  down”  by  a  “candidate  location 
flag.”  Figure  5.3  shows  the  received  signal  ratios  (constant  level  line)  and  the 
corresponding  candidate  location  flags,  for  four  antenna  pairs. 

Again,  we  should  remember  that  the  real  comparison  is  being  made 
between  many  more  points,  since  there  is  a  third  dimension  (height).  For  simplicity, 
we  show  here  only  one  of  36  height  slices. 
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The  dots  in  Figure  5.3  (upper  portion)  represent  the  real  calculated 
values,  while  the  connecting  curves  axe  made  by  extrapolation.  Therefore,  a  can¬ 
didate  location  flag  is  set  only  when  the  measured  signal  ratio  is  within  a  specified 
distance  from  the  calculated  signal  ratio.  This  distance  is  determined  by  the  width 
of  a  comparison  window  which  is  a  software  parameter.  The  wider  the  window,  the 
more  candidate  location  flags  we  will  set.  On  the  other  hand,  we  cannot  allow  too 
narrow  a  window,  since  there  is  a  risk  of  loosing  correct  locations.  It  is  obvious  that 
as  we  use  higher  resolution  in  the  lookup  table,  we  can  also  narrow  the  comparison 
window. 

Figure  5.4  shows  the  candidate  location  flags  for  the  same  scenario  as 
in  Figure  5.3,  but  with  a  wider  comparison  window.  Notice  that  more  flags  appear 
in  this  case. 

For  a  higher  sea  state,  there  are  random  fluctuations  in  signal  level, 
and  the  comparison  becomes  harder.  We  must  open  the  window  wider,  but  at  the 
same  time  more  incorrect  locations  are  flagged.  A  practical  compromise  is  required. 

f.  Validation: 

After  we  have  all  of  the  candidate  location  flags,  we  need  to  decide 
which  location  is  correct.  The  method  of  doing  this  is  by  doing  a  cross  location: 
comparing  between  all  the  flags,  and  checking  which  one  repeats  in  all  of  the  ratios. 

If  after  this  process  there  is  more  than  a  single  location  (more  than 
only  one  cross  location),  the  measurement  cannot  give  a  correct  location,  so  we 
discard  it,  and  start  a  new  session. 

Practically,  using  three  receiving  antennas  (two  ratios)  gives  no  single 
location  in  most  cases,  and  therefore  we  need  to  add  more  antennas  in  order  to  have 
more  ratios.  During  the  experiments  we  tried  using  different  numbers  of  antennas, 
and  the  final  sufficient  result  uses  five  receiving  antennas. 
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Figure  5.5  shows  the  candidate  location  flags  for  each  antenna  pair  and 
the  crossing  between  them,  giving  a  single  valid  location  marked  by  a  diamond.  The 
sea  state  is  zero,  and  the  window  width  is  2.5%.  Here,  we  can  see  not  only  a  single 
“slice”,  but  the  whole  combination  of  range  and  height.  From  this  it  is  obvious  that 
a  single  pair  of  antennas  gives  no  information  at  all  (too  many  candidate  location 
flags). 

For  the  case  of  three  antennas,  (see  Figure  5.6),  the  ambiguity  is  re¬ 
duced  significantly  (only  two  flags  overlapping),  while  adding  a  fourth  antenna, 
reduces  the  location  to  a  single,  correct  one.  In  this  particular  case,  four  receiving 
antennas  give  the  correct  location.  In  many  cases,  though,  four  antennas  will  still 
leave  some  order  of  ambiguity,  but  the  fifth  one  will,  in  most  cases,  give  a  correct, 
single  location. 

The  ambiguity  level  reduction  as  a  function  of  number  of  receiving 
antennas  is  well  illustrated  in  Figure  5.6,  for  sea  state  zero. 

For  a  wider  window  there  will  be  more  flags,  as  shown  in  Figure  5.7, 
and  here  we  can  see  that  only  the  fifth  antenna  solves  the  ambiguity. 

When  the  sea  state  is  higher,  a  wide  window  will  lead  to  a  large  level 
of  ambiguity.  Figure  5.8  shows  a  12%  window  in  sea  state  3,  where  there  are  four 
locations.  Reducing  the  width  to  5%  will  reduce  the  number  of  locations  to  two,  as 
shown  in  Figure  5.9,  while  the  ambiguity  will  totally  disappear  for  a  2.5%  window, 
as  in  Figure  5.10. 

g.  Output: 

The  previous  step  brought  up  a  single  location  flag,  which  shows  the 
correct  location  of  the  transmitter  in  the  range-height  table.  The  output  step 
translates  the  table  location  into  the  real  coordinates  of  range  and  height,  and 
types  it  as  an  output  result. 
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C.  FINAL  RESULTS 

As  mentioned  earlier  in  this  chapter,  the  fined  resolution  used  here,  is  130  steps 
of  50  meters  each,  from  the  minimum  range  up  to  7.5  Km,  and  then  125  steps  of 
100  meters  each  to  the  maximum  range.  The  comparison  window  width  is  2.5%. 

With  these  parameters,  estimation  results  are  given  in  Figures  5.11  and  5.12, 
for  sea  states  0  and  3,  respectively.  It  should  be  mentioned  that  the  estimated 
range  sometimes  deviates  from  the  actual  range  by  a  few  tens  of  meters,  which 
cannot  be  seen  in  the  graph,  due  to  resolution. 
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Figure  5.3:  Comparison  Process  with  2.5%  Window. 
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Figure  5.5:  Validation 
Diamond. 


CANDIDATE  LOCATION  FLAGS  FOR  2.5  PERCENT  WINDOW 

2  RECEIVING  ANTENNAS  3  RECEIVING  ANTENNAS 
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Figure  5.6:  Ambiguity  Level  vs.  Number  of  Receiving  Antennas;  Sea 
State  0  and  2.5%  Window. 
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VALIDATED  LOCATION  FLAGS  FOR  12  PERCENT  WINDOW 
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Figure  5.7:  Ambiguity  Level  vs.  Number  of  Receiving  Antennas;  Sea 
State  0  and  12%  Window. 
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Figure  5.8:  Ambiguity  Level  vs.  Number  of  Receiving  Antennas;  Sea 
State  3  and  12%  Window. 
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CANDIDATE  LOCATION  FLAGS  FOR  5  PERCENT  WINDOW 

2  RECEIVING  ANTENNAS  3  RECEIVING  ANTENNAS 
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Figure  5.9:  Ambiguity  Level  vs.  Number  of  Receiving  Antennas;  Sea 
State  3  and  5%  Window. 
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VALIDATED  LOCATION  FLAGS  FOR  2.5  PERCENT  WINDOW 

2  RECEIVING  ANTENNAS  3  RECEIVING  ANTENNAS 
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Figure  5.10:  Ambiguity  Level  vs.  Number  of  Receiving  Antennas;  Sea 
State  3  and  2.5%  Window. 
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ESTIMATED  VERSUS  ACTUAL  RANGE  -  SEA  STATE  0 


(&313IAI)  3  0  N  V  y 


Figure  5.11:  Estimated  Results  (Range);  Sea  State  0. 
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Figure  5.11  (continued):  Estimated  Results  (Height);  Sea  State  0. 
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Figure  5.12:  Estimated  Results  (Range);  Sea  State  3. 
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D  VERSUS  ACTUAL  HEIGHT  -  SEA  STATE 


VI.  ADDITIONAL  PARAMETERS  CONSIDERATION 


This  chapter  deals  with  some  parameters  that  have  influence  on  the  calcula¬ 
tions  and  terms  used  for  the  estimation  process.  A  few  of  them  were  mentioned 
but  not  discussed  for  ail  cases,  and  the  others  were  not  discussed  previously. 

A.  VERTICAL  POLARIZATION 

So  far  we  dealt  with  horizontal  polarization,  and  all  the  calculations  were  done 
for  this  polarization.  Vertical  polarization  complicates  the  problem  a  bit,  since  the 
reflection  coefficient,  p,  changes  significantly  with  the  grazing  angle,  as  shown  in 
Figure  3.4. 

Let  us  check  some  representative  cases  in  order  to  see  how  much  the  reflection 
coefficient  really  changes  due  to  vertical  polarization. 

Case  1:  transmitter  is  located  10  Km  apart  from  the  receiver,  at  a  height  of 
10  m,  and  it  transmits  with  a  frequency  of  10  GHz.  The  receiver  antennas  are  at 
heights  of  12  m  and  20  m.  So,  our  case  has  the  parameters: 

R  =  10  km, 
ht  =  10  m, 

/  =  10  GHz  (A  =  3  cm), 
h rl  -  12  m,  and, 
hr2  =  20  m. 

Using  equations  (5.4)  to  (5.12),  we  find  that  the  grazing  angles  are  equal  to: 
il'i  —  0.126°,  and  xl'-i  =  0.172°,  respectively  for  the  receiving  antenna  heights. 
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Using  the  correct  graph,  (Figure  3.4),  we  find  that  the  reflection  coefficient  due  to 
the  reflection,  equals  0.98  for  both  grazing  angles. 

Case  2:  in  the  second  case  the  transmitter  is  higher,  and  the  range  between 
the  receiver  to  the  transmitter  is  reduced: 

R  =  4  km, 
ht  =  40  m, 

/  =  10  GHz  (A  =  3  cm), 
hr\  =  12  m,  and, 
hT 2  =  20  m. 

Using  the  same  method  as  in  the  first  case,  we  find  here,  that  the  reflection 
coefficient  for  the  lower  receiving  antenna  is  0.82,  while  for  the  higher  antenna  it 
is  0.80. 

Case  3:  here  we  take  an  extreme  case  of  very  high  transmitter,  as  well  as 
receiver  antennas: 

R  =  4  km, 
ht  =  100  m, 

/  =  10  GHz  (A  =  3  cm), 
hri  =  38  m,  and, 
hr2  =  50  m. 

In  this  case,  the  reflection  coefficients  are  0.58  for  the  38  m  receiving  antenna, 
and  0.55  for  the  50  m  antenna. 

We  can  conclude  that  for  low  geometry  (transmitter  and  receiver),  and  large 
ranges,  the  reflection  coefficient  can  be  assumed  the  same  for  all  receiving  antennas. 
For  shorter  ranges  and  larger  heights,  we  should  take  into  account  the  changes  in 
the  reflection  coefficient  for  the  different  receiving  antennas. 
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Practically,  for  the  vertical  polarization  case,  we  need  to  modify  the  off  line 
tables  by  the  amount  of  the  reflection  coefficient  changes.  This  can  be  done,  since 
we  know  how  to  calculate  it  with  the  given  parameters,  and  for  the  whole  range 
and  transmitter  height  field.  The  only  parameter  that  has  doubt  is  the  kind  of 
polarization.  This  verification  must  be  done  by  another  system,  and  provided  as 
an  input  to  the  estimation  process. 

B.  RECEIVER  HEIGHT 

As  we  saw  in  Chapter  4,  the  higher  the  receiver  antenna  is,  the  more  complex 
the  curve  is.  In  our  work  we  used  a  maximum  receiving  antenna  height  of  20  meters, 
since  it  is  a  reasonable  height  for  a  conventional  missile  boat.  For  a  larger  ship  (such 
as  a  corvette  or  destroyer,  for  example),  we  shall  need  higher  receiving  antennas. 
The  only  influence,  as  mentioned  before,  is  on  the  curves  complexity.  The  solution 
for  this  is  to  use  a  greater  resolution  in  the  look  up  tables,  to  compensate  for  the 
higher  ambiguity. 

C.  SHIP  TILTING 

For  the  ships  of  our  interest  (relatively  small  missile  boats),  the  worst  case 
rolling  angle  can  be  around  15°  for  cross  wind  cruising  direction  in  sea  state  3. 
For  sea  state  4,  the  extreme  angle  is  about  20°  for  the  same  conditions.  Since 
the  tilt  occurs  for  both  sides,  we  can  say  that,  to  a  first  degree  of  approximation, 
the  average  tilt  is  0°.  Even  though,  let  us  check  the  worst  case.  The  geometry  of 
the  problem  is  given  in  Figure  6.1.  We  will  now  check  the  tilt  impact  on  the  look 
up  tables  for  two  cases:  the  closest  antennas  (12  meters  and  14  meters),  and  the 
extreme  antennas  (12  meters  and  20  meters). 
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Figure  6.1:  Tilted  Ship  Geometry. 

Referring  to  Figure  6.1,  in  our  case,  h\  =  12  m,  h2  =  14  rp  and  h$  =  20  m. 
The  maximum  tilt  (roll)  angle,  amAX,  is  taken  as  15°  for  sea  state  3,  and  20°  for 
sea  state  4. 

From  Figure  6.1  we  see  that  the  actual  height  of  each  receiving  antenna,  after 
being  tilted,  is: 

h\  =  h{  ■  cos  a  ■  (i  =  1,2,  •  •  •  ,5).  (6.1) 

For  sea  state  3  (amax  =  15°),  we  obtain: 

h\  =  12  -cos  15°  =  11.6  m  (where  h\  =  12  m),  (6.2) 

and 

h'5  =  20  ■  cos  15°  =  19.3  m  (where  h2  =  20  m).  (6.3) 

The  actual  distance  between  the  antennas  is  7.7  meters,  instead  of  a  theoretical 
value  of  8  m.  Figure  6.2  shows  the  ratio  Tjj/rjj  for  receiving  antennas  located 
perpendicular  to  the  sea  level,  in  heights  of  20  meter  and  12  meter,  respectively. 
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Figure  6.3:  Ratio  Between  Ratios  of  Attenuation  Coefficients;  Sea  State 
3,  Receiving  Antennas  at  12  m  and  20  m. 


Now  we  will  repeat  the  calculations  for  receiving  antennas  at  heights  of  14 
meters  and  12  meters.  Note  that  h\  remains  the  same  (11.6  m),  while  h'2  becomes: 


h'2  =  14  •  cos  15°  =  13.5  m  (where  /12  =  14  m).  (6.4) 
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Figure  6.4:  Ratio  Between  Ratios  of  Attenuation  Coefficients;  Sea  State 
3,  Receiving  Antennas  at  12  m  and  14  m. 

The  actual  distance  here  is  1.9  meters,  instead  of  2  meters.  Figure  6.4  shows 
the  ratio  between  rjj /r/j  (perpendicular  antennas),  and  (»7| ) VI7??  (tilted  anten¬ 
nas). 

Repeating  the  same  calculations  for  sea  state  4  (amax  =  20°),  we  obtain: 

h[  =  12  •  cos 20°  =  11.3  m  (where  hi  =  12  m),  (6.5) 

and 

h'5  =  20  ■  cos 20°  =  18.8  m  (where  h5  =  20  m).  (6.6) 
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The  actual  distance  between  the  antennas  is  7.5  meters,  instead  of  a  theoretical 
value  of  8  meters.  Figure  6.5  shows  both  ratios  of  r]l/rj\  and  (775 )2 / (77^  )2 ,  while 
Figure  6.6  shows  the  ratio  between  those  two  ratios. 

For  the  closest  antennas  we  obtain: 

h'2  =  14  •  cos  20°  =  13.2  m  (where  h 2  =  14  m),  (6-7) 

while  h[  remains  the  same,  11.3  m.  The  actual  difference  is,  therefore,  1.9  meters, 
instead  of  2  meters.  Figure  6.7  shows  the  ratio  between  (perpendicular 

antennas),  and  (v'2)2 /(Vi)2  (tilted  antennas). 

To  summarize  these  results,  we  see  that  the  method  is  sensitive  to  ship  tilting, 
since  there  is  a  difference  between  the  tables  for  perpendicular  and  tilted  config¬ 
urations.  The  conclusion  of  these  results  is  that  the  receiving  antennas  must  be 
stabilized. 

D.  RECEIVER  NOISE 

In  the  case  of  our  work,  the  received  signal  to  noise  ratio  was  assumed  to  be 
relatively  high.  The  ranges  of  interest  are  less  than  20  Km,  while  the  missile  ERP 
is  in  the  order  of  100  to  120  dBm.  In  addition,  the  transmitted  signal  propagates 
only  one  way.  Therefore,  any  received  signal  is  much  larger  than  the  receiver  noise, 
and  the  signal  to  noise  ratio  is  high  enough,  so  we  can  disregard  the  noise. 

E.  ATMOSPHERIC  ATTENUATION 

The  atmospheric  attenuation  can  be  treated  in  two  ways:  modify  the  look  up 
tables  with  the  attenuation  data,  since  all  the  relevant  parameters  are  known,  or. 
more  practically,  to  disregard  it.  Since  we  deal  with  ratios,  the  attenuation  will 
have  an  influence  only  if  the  path  difference  is  significant.  As  we  already  discussed 
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Figure  6.5:  Ratio  Between  Attenuation  Coefficients;  Sea  State  4, 
Receiving  Antennas  at  12  m  and  20  m. 
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Figure  6.6:  Ratio  Between  Ratios  of  Attenuation  Coefficients;  Sea  State 
4,  Receiving  Antennas  at  12  m  and  20  m. 

in  Chapter  III,  the  path  difference  is  negligible  compared  to  the  total  path  length. 
Therefore,  the  attenuation  is  practically  the  same  for  the  direct  and  indirect  path, 
so  it  has  no  influence  on  the  problem,  hence  can  be  neglected. 


F.  OTHER  MODELS 

Appendix  A  gives  other  models  for  reflection  coefficients,  instead  of  the  Ament 
model,  that  we  used  in  this  work.  To  use  another  model,  the  only  thing  that  should 
be  taken  care  of  is  modifying  the  model  equation  in  the  computer  program.  This 
point  will  be  discussed  in  detail  in  the  next  chapter. 


75 


Figure  6.7:  Ratio  Between  Attenuation  Coefficients;  Sea  State  4, 
Receiving  Antennas  at  12  m  and  14  m. 

G.  EARTH  CURVATURE 

In  our  work,  we  assumed  a  plane  earth.  In  general,  some  geometries  need 
to  take  account  of  the  earth’s  curvature.  Our  particular  geometry  results  in  no 
influence  from  the  earth  curvature,  as  shown  below.  The  extreme  condition  is  for 
the  lowest  receiving  antenna,  and  the  lowest  missile. 

The  maximum  range  that  still  allows  plane  earth  assumption  is  calculated  by 
using  equation  6.8  [Ref.  7].  The  problem  geometry  is  shown  in  Figure  6.8. 

Rma x  =  \j2aehr  +  hf.  +  yj 2 aeht  4-  h\ ,  (6.8) 
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Figure  6.8:  Plane  Earth  Assumption  Geometry. 

where: 

ae  is  the  effective  earth  radius  (=  4/3  x  Re), 

Re  is  the  earth  radius  (=  6370  Km), 
hr  is  the  receiving  antenna  height,  and, 
ht  is  the  transmitter  height. 

For  all  practical  purposes,  the  term  of  h 2  under  the  square  root  can  be  ne¬ 
glected. 

Our  extreme  case  is  when  the  missile  cruises  in  ht  =  5  m,  and  the  lowest 
receiving  antenna  is  located  at  height  of  ht  =  12  m.  For  this  geometry  we  obtain 
a  maximum  range  of: 

RmiX  =  V2  ■  8493  ■  103  •  12  +  V2  ■  8493  ■  103  •  5  =  23.5  Km. 
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Since  we  deal  with  relative  ranges  of  up  to  20  Km,  we  can  disregard  the  earth 
curvature. 
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VII.  COMPUTER  PROGRAM 


The  computer  program  in  this  work  was  written  in  APL.  The  listing  is  given  in 
Appendix  B  and  it  includes  all  the  necessary  comments.  This  chapter  presents  an 
overview  of  the  whole  program,  and  then  explains  each  function  in  greater  detail. 

The  program  consists  of  a  few  functions  which  call  each  other  in  sequence,  and 
a  main  management  function. 

The  sequence  follows  exactly  after  the  one  described  in  Chapter  V,  under 
“Detailed  Discussion.”  All  the  inputs  that  must  be  given  by  external  means,  are 
inserted  here  by  the  user,  &s  software  parameters. 

The  program  is  started  by  typing  RUN.  The  user  is  asked  to  enter  sea  state 
(0  to  7),  and  with  this  data,  the  RUN  function  calls  the  main  function  for  receiving 
data,  named  DATAIN.  This  function  automatically  creates  more  needed  parame¬ 
ters,  in  a  random  manner,  between  given  limits.  It  is  also  possible  for  the  user  to 
dictate  each  one  of  these  parameters,  instead  of  creating  them  randomly  by  the 
program: 

•  frequency  -  between  2  to  18  GHz, 

•  ERP  (Effective  Radiated  Power)  -  between  100  to  120  dBm  (107  to  109 
Watt), 

•  initial  range  for  simulation  (range  where  missile  starts  to  radiate)  -  between 
10  to  20  Kilometers,  and, 

•  transmitter  height  -  between  5  to  40  meters. 

Other  relevant  parameters  are  programmed  off-line  as  software  parameters, 
and  can  be  changed  easily  at  any  time.  Those  are: 


•  window  width,  as  discussed  in  Chapter  V, 

•  minimum  transmitter  height  for  the  look  up  table, 

•  number  of  steps  for  the  transmitter  height  for  the  table, 

•  resolution  of  transmitter  height  steps  for  the  table, 

•  minimum  range  for  simulation  (if  the  missile  passed  this  range  without  being 
intercepted,  the  simulation  is  failed  and  stops), 

•  range  resolution  for  the  look  up  table,  both  high  resolution  for  short  range, 
and  regular  resolution  for  longer  range, 

•  number  of  range  steps  for  each  range  resolution, 

•  resolution  of  simulation  steps  (range  distance  between  each  receiving),  and, 

•  locations  (heights)  of  receiving  antennas. 

The  function  then  calculates  the  maximum  range,  and  shows  the  major  pa¬ 
rameters  on  the  screen  (a  full  session  is  demonstrated  later  in  this  chapter.) 

With  all  the  necessary  data,  DATAIN  returns  to  the  management  function, 
RUN.  In  the  next  step,  RUN  sends  all  the  relevant  parameters  to  the  main  function 
STARTPRO,  and  activates  it. 

STARTPRO  is  responsible  for  all  the  major  calculations: 

•  creating  off  line  tables  which  include  also  the  attenuation  coefficient,  p2 , 
and  the  reflection  coefficient,  p, 

•  receiving  five  signals  in  the  receiving  antennas,  and  divide  their  powers  to 
determine  the  ratios, 

•  comparing  between  the  off  line  tables  and  the  actual  ratios, 

•  marking  the  correct  (single)  location,  and  finally, 

•  sending  the  location  to  the  management  function. 
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The  signals  for  STARTPRO  are  created  by  the  simulator  SIM.  This  function 
is  activated  by  STARTPRO,  every  time  it  needs  a  new  set  of  receiving  signals. 

SIM  receives  the  relevant  parameters  from  DATAIN,  and  creates  five  signals 
due  to  those  parameters  and  the  geometry.  The  signal  is  calculated  using  equa¬ 
tion  7.1. 

S,  =  ERP-e-,^,  (7.1) 

where: 

ERP  is  the  effective  radiated  power  of  the  transmitter, 

R  is  the  relative  range  between  the  transmitter  and  the  receiver, 

t}2  is  the  attenuation  coefficient,  calculated  by  using  Equation  (3.5),  and, 

p  is  the  reflection  coefficient,  calculated  by  using  Equation  (3.20). 

After  the  management  function  RUN  received  the  location  of  the  transmitter 
from  STARTPRO,  it  calls  the  translating  function  RHFIND  in  order  to  translate 
the  coordinates  of  the  table  into  range  and  height,  and  to  type  it. 

A  conceptual  flow  chart  is  given  in  Figure  7.1,  and  a  simulation  session  is 
shown  in  Figure  7.2. 
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Figure  7.1:  Estimation  Process;  Conceptual  Flow  Chart. 


82 


♦RUN 

ENTER  SEA  STATE: 


Figure  7.2: 


0  -  CALM  1  -  SMOOTH  2  -  SLIGHT 

3  -  MODERATE  4  -  ROUGH  5  -  VERY  ROUGH 

6  -  HIGH  7  -  VERY  HIGH 


♦2  (for  example) 

RANGE  INTERVAL  FOR  SIMULATION: 

MINIMUM  RANGE  (IN  METERS):  1000 
MAXIMUM  RANGE(IN  METERS):  20000 
SEA  STATE: 

2 

FREQUENCY  (IN  GHZ): 

8.85 

INITIAL  RANGE  (IN  KILOMETERS): 

14.5 

TRANSMITTER  HEIGHT  (IN  METERS)  -  ACTUAL: 
36 

RECEIVER  ANTENNAS  HEIGHTS  (IN  METERS): 

12  14  16  18  20 

TRANSMITTER  RANGE  (IN  KILOMETERS): 

14.4 

TRANSMITTER  HEIGHT  (IN  METERS): 

36 

UPDATED  ACTUAL  RANGE  (IN  KILOMETERS): 
14.4 


Simulation  Session  Example.  (Lines  Starting  with  Asterisk 
(*)  are  Entered  by  the  User.) 
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Vm.  SUMMARY  AND  CONCLUSIONS 


The  work  which  had  been  done  during  this  thesis  research,  together  with  the 
results  obtained,  show  that  the  suggested  method  for  passive  range  estimation 
using  multipath  effect  can  probably  be  used.  The  work  developed  the  method 
theoretically,  and  proved  its  feasibility. 

The  accuracy  of  the  estimated  range  and  height  depends  upon  many  param¬ 
eters,  as  well  as  receiving  antennas  configuration.  A  trade  off  should  be  done  be¬ 
tween  number  of  receiving  antennas,  their  heights,  separation  and  window  width, 
as  described  in  Chapter  V. 

There  is  no  analytical  way  for  finding  the  optimum  trade  off,  and  only  trial 
and  error  sessions  can  achieve  a  local  optimum,  as  was  done  during  this  research. 

Another  question  of  interest  is  the  accuracy  of  the  rough  sea  model  (Ament 
model  in  our  case). 

The  answer  to  this  question  can  only  be  obtained  by  sea  trials.  Those  trials 
should  answer  the  questions  listed  below: 

•  Is  the  theoretical  analysis  of  the  additional  parameters  in  Chapter  VI  cor¬ 
rect,  or  should  we  modify  these  parameters? 

•  What  is  the  estimation  accuracy  if  an  incorrect  sea  state  is  entered,  and 
is  there  a  preferred  direction  for  entering  sea  state  (is  it  better  to  enter  a 
higher  sea  state  in  case  of  uncertainty,  or  lower  state)? 

•  How  does  the  wind  direction  influence  the  accuracy,  if  at  all,  (upwind,  down¬ 
wind,  and  crosswind)? 

•  What  is  the  practical  feasibility  of  the  suggested  method,  and  with  which 
limitations  (such  as  sea  state  and  wind  direction)? 
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We  should  notice  that  the  Ament  model  neglects  the  shadowing  effect,  al¬ 
though  it  deals  with  small  grazing  angles,  where  this  effect  does  occur.  It  is  ques¬ 
tionable  how  far  the  shadowing  effect  changes  the  reality  from  the  model. 

A  possible  influence  of  the  shadowing  effect  is  the  fact  that  the  amplitude  of 
the  reflected  wave  will  be  lower,  so  that  the  signal  minima  will  not  be  so  sharp  as 
is  shown  in  Figure  4.2.  A  decrease  in  tb  amplitude  of  the  specular  reflection  will 
also  occur  occur  because  of  high  sea  state,  and  capillary  waves.  If  sea  trials  show 
consistent  errors,  we  can  conclude  that  the  Ament  model  is  inadequate.  In  this 
case  it  will  be  necessary  to  develop  another  more  accurate  model. 

It  can  be  predicted  that  for  sea  state  0,  the  results  will  be  more  accurate  than 
in  any  other  case,  while  the  more  roughness,  the  less  accurate  the  results. 

In  order  to  improve  the  estimation  process,  and  to  reach  a  lower  rate  of  errors, 
we  can  modify  the  program.  The  modified  program  will  consider  also  the  past 
readings,  and  compare  each  present  estimate  with  the  previous  ones.  If  the  range 
becomes  shorter  by  a  reasonable  increment,  and  the  transmitter  height  remains 
almost  the  same,  we  can  conclude  that  the  estimate  is  correct.  On  the  other  hand, 
if  the  estimate  does  not  meet  this  test,  then  another  estimate  should  be  made,  and 
compared,  and  so  on. 
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APPENDIX  A:  REFLECTION  COEFFICIENT  MODELS 


In  Chapter  III  we  discussed  the  most  popular  model  for  reflection  coefficient 
of  a  rough  sea,  the  Ament  model.  There  axe  several  other  models,  which  are  given 
below.  It  should  be  noted  that  there  are  additional  models  which  we  will  not  cover 
here. 


A.  BULLINGTON  MODEL 


This  model  distinguishes  between  different  surface  profiles  [Ref.  3,12]: 

1.  Sawtooth  Profile 

The  wave  height  is  distributed  uniformly,  but  does  not  necessarily  have 
uniform  slope.  The  phase  values  have  equal  probability  between  ~(pm  and  +<pm , 
where 


4* m  — 


27r<7fc  sin  ip 
A  : 


(A.l) 


a/,  is  the  RMS  deviation  of  the  waves  height, 
ip  is  the  grazing  angle,  and, 

A  is  the  signal  wavelength. 


The  reflection  coefficient  for  the  sawtooth  profile  is  given  by: 


p  —  _____  —  smc0m, 

<Pm 


and  is  shown  in  Figure  A.l. 


(A.2) 
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Figure  A.l:  Bullington  Model  -  Sawtooth  Profile. 


2.  Sinusoidal  Profile 

The  reflection  coefficient  in  this  case  is  given  by  Equation  (A. 3),  and  is 
shown  in  Figure  A. 2. 


(  27T (t h  sin 


where  Jo  is  the  Bessel  function  of  the  first  kind  and  zero  order. 
3.  Triangular  Profile: 

Here,  the  reflection  coefficient  is  given  by: 

2 


__  /  sin  <f>i  \ 

"  W*  )  ’ 


(A.  3) 


(-4.4) 
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Figure  A.2:  Bullington  Model  -  Sinusoidal  Profile. 

where  the  maximum  phase  is: 


4>t  = 


n(7h  sin  0 
A 


(.4.5) 


Figures  A. 3  and  A. 4  show  the  triangular  profile  reflection  coefficient. 

4.  Sine  Square  Profile 

The  reflection  coefficient  is  more  complex  in  this  model,  and  is  given  by 
Equation  (A.6),  and  shown  in  Figure  A. 5. 


p  =  cos 4>i  ■  J0(<f>t), 

where  <t>t  is  the  same  as  given  in  Equation  (A. 5). 


(.4.6) 
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SIGMA  x  XI  4-  LAMBDA  (MILRAD) 


Figure  A.3:  Bullington  Model  -  Triangular  Profile. 

B.  DAVIES  MODEL 

The  Davies  model  assumes  that  the  surface  has  a  random  behavior,  which 
yields  a  gaussian  height  distribution  and  correlation  coefficient.  This  model  neglects 
the  shadowing  effect,  and  is  good  for  small  to  medium  grazing  angles  (ip  <  40°).  In 
order  to  describe  the  surface,  we  introduce  an  arbitrary  parameter,  m/a.  Experi¬ 
ments  above  the  sea  surface  show  that  a  value  of  m/a  =  50  to  60,  gives  sufficient 
correlation  between  theoretical  and  experimental  results.  [Ref.  3]. 

The  reflection  coefficient  for  this  model  is  given  by: 


4ec  ip  ■  exp  tan2  ip  . 


(A- 7) 


Figure  A. 6  shows  the  reflection  coefficient.  The  value  m/2a  is  27  (m/a  =  54). 


89 


Figure  A. 4:  Bullington  Model  -  Triangular  Profile  (expanded). 
C.  ISAKOVICH  MODEL 


In  this  model  we  assume  that  the  surface  is  irregular,  isotropic  in  all  directions, 
and  has  a  normal  height  distribution.  The  parameter  m/a  describes  the  surface  (as 
in  the  Davies  model).  It  is  the  only  parameter  which  dictates  the  intensity  of  the 
scattering  in  a  given  direction.  Changing  the  ratio  will  change  the  profile  shape, 
as  shown  in  Figure  A. 7  (m/a  =  10),  and  Figure  A. 8  ( m/a  =  20).  The  shadowing 
effect  is  neglected  and  only  fear-field  is  calculated.  The  Isakovich  model  calculates 
the  dispersion  of  waves  for  random  surface  [Ref.  3]: 


m 

P  =  — 
a 


1  —  cos  ip  cos  Et  cos  9  +  sin  xp  sin  Et 


(simp  -I-  sin  Et)2 


■  exp  — 


1  /  m  \  2  cos2xp  +  cos2  Et  —  2  cos  ip  cos  Et  cos  9 
8  \  a  )  (simp  +  sinF,)2 


(-4.8) 
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Figure  A. 5:  Bullington  Model  -  Sine  Square  Profile. 

This  model  is  much  more  complicated  than  the  others,  since  it  has  two  more 
parameters: 

Et  is  the  observation  angle,  and, 

9  is  the  azimuth  angle. 


In  order  to  illustrate  results  comparable  to  those  for  the  other  modes,  we 
assume  that  the  azimuth  angle  is  zero  (9  =  0).  For  simplicity,  we  calculate  the 
reflection  coefficient  for  two  cases:  zero  observation  angle,  Et  =  0°  (transmitter 
and  receiver  at  the  same  level),  and  Et  =  —30°  (transmitter  higher  than  the 
receiver). 
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GRAZING  ANGLE 


(MILRAD) 


Figure  A. 6:  Davies  Model. 


Using  the  assumption  9  =  0,  we  simplify  Equation  (A. 8)  to: 


This  case  is  illustrated  in  Figures  A. 7  and  A. 8  for  two  different  values  of 
m/a.  For  —30°  observation  angle  Equation  (A. 11)  applies.  The  result  is  shown  in 
Figure  A. 9. 


0  200  400  600  800  1000 


GRAZING  ANGLE  (MILRAD) 


Figure  A. 7:  Isakovich  Model  with  Zero  Observation  Angle  ( Et  =  0°), 
mja  =  10. 


Figure  A.8:  Isakovich  Model  with  Zero  Observation  Angle,  ( Et  —  0°), 
mja  —  20. 


t 
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REFLECTION  COEFFICIENT 


O 


Figure  A. 9:  Isakovich  Model  with  Observation  Angle  of  —30°,  m/cr  =  20 
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APPENDIX  B:  COMPUTER  PROGRAM  LISTING 


This  appendix  gives  the  listing  of  the  computer  program  described  in  Chapter 
VII.  RUN  is  the  management  function,  which  calls  the  other  functions  in  order. 
DATAIN  is  the  function  responsible  for  the  input  data  for  further  calculations. 
STARTPRO  is  the  main  function  which  does  all  the  major  calculations,  as  well  as 
creating  the  look-up  tables,  and  comparing  the  ratios  of  the  actual  received  signals 
with  the  calculated  ratios.  SIM  is  the  simulation  function,  which  simulates  the 
received  signals  from  the  missile.  RHFIND  is  the  function  which  translates  the 
location  of  the  missile  in  the  table  into  a  range  and  height. 


**RUN 


[0] 

RUN 

LI] 

a  MANAGMENT  FUNCTION 

[2] 

'ENTER  SEA  STATE-.' 

[3  ] 

i  i 

[4] 

' 0-CALM  1 -SMOOTH  2-SLICHT' 

[5] 

•3 -MODERATE  4 -ROUGE  5 -VERT 

ROUGH 

C6] 

'6-RICH  7 -VERY  HICH' 

[7] 

■  i 

L8] 

SEA-D 

[9] 

a  ENTERING  ROUCHNESS  PARAMETER 

CIO] 

i  i 

[11] 

DATAIN 

C12] 

a  INITIALIZATION 

[13] 

RR-.F  STARTPRO  HR3 

C 1 4  ] 

a  STARTINC 

-C15] 

*(FLAG= 1 )/ST0 

[16] 

a  END  OF  SESSION  IF  RANCE  CANNOT  BE 

FOUND 

[17] 

RHFIND 

[18] 

a  CO  TO  CHECK  REAL  RANCE  AND  HEIGHT 

-C19] 

STO : -0 

**DATAIN 

[0]  DATAIN 

Cl]  a  INITIALIZATION  FUNCTION 

[2]  r  aooaon oo 

[ 3 ]  MINHEIGT** 

[4]  r  1 MINHEIGT '  IS  THE  MINIMUM  TX.  HEIGHT  (.LESS  1  METER) 

[5]  r  00000000 

[63  HEICTST+Z6 

[73  p  'HEIGTST'  IS  THE  NUMBER  OF  HEIGHT  STEPS 

[8]  r  aaooaooo 

[9]  HEIGTRES* 1 

[103  R  ' HEIGTRES '  IS  THE  RESOLUTION  OF  HEIGHT  STEPS 

[113  a  ooooooaa 

[123  00.975 

[133  a  aoaooaoa 

Cl4]  01.025 

[153  R  1 C 1 , ' D'  ARE  LONER  AND  UPPER  LIMITS  FOR  OK  (START PRO) 

[16]  a  ODDODODQ 

[17]  RMIN* 1000 

[183  fl  ' RMIN ’  IS  MINIMUM  RANGE 

[193  RMINS+RMIN 

C203  a  'RMIN S'  IS  MINIMUM  RANGE  FOR  SIMULATION 

[2i3  r  aaooaooo 

C22]  N STEPS* 125 

[233  r  'NSTEPS'  IS  THE  NUMBER  OF  RANGE  STEPS 

[24]  a  aaoooooa 

[25]  RESO* 100 

C26]  fl  'RESO'  IS  TRE  RANCE  STEPS  RESOLUTION  (LENGTH  OF  EACH  STEP ) 

[27]  o  oooooaaa 

C2S3  fiNSTEPS'*’  130 

[29]  fl  'HN STEPS'  IS  THE  NUMBER  OF  RICH  RESOLUTION  RANGE  STEPS 

[30]  fl  OOODODDO 

[31]  HRESO*50 

[32]  fl  'HRESO'  IS  THE  RANCE  STEPS  (HIGH)  RESOLUTION  FOR  SMALL 

[33]  A  RANGES 

[34]  «  aoooaoaa 

[35]  SIMSTEP* 50 

[36]  a  'SIMSTEP'  IS  THE  RANGE  STEPS  RESOLUTION  IN  THE  SIMULATION 

[3  7]  a  FUNCTION  (i7 SIM) 

[383  F* (1999+(?17001))*1000 

[393  a  oaaaaaoa 

[403  A09.35 

[41]  a  <F<  IS  FREQUENCI  (IN  GHZ);  2*18  GHZ 
[423  £7?010000000x?100 

[433  a  aaoooooa 

[44]  aERP*100000000 

[45]  a  'ERP'  IS  TRANSMITTER  EFFECTIVE  RADIATED  POWER  (IN  WATTS); 

[46]  a  1E7+1E9  WATT  (100+120  DBM) 

[47]  RI+2QQQ0* (500  +  17300  )  + 1000 

[48]  a  'RI'  IS  INITIAL  RANCE  (IN  METERS);  10,000*16,000  METER 

[49]  a  OOOOOOOQ 
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□□CDDODO 


[50]  o  RI  -*-15400 

[51]  HAT*-1:  +  ?  3  6 

[52]  a 

[53]  aHAT+22 

[5u]  a  'HAT'  IS  TRANSMITTER  HEIGHT  (.IN  METERS );  5  +  uo  METER 

[55]  o,V/?3*-l  4  +  5  ?  16 

[56]  a  'HR 3'  ARE  RECEIVER'S  ANTENNAS  HEIGHTS  (IN  METERS) 

[57]  a  'HR 3'  IS  BETWEEN  15+30  METERS 

c  5  s ]  a  aaaaoooa 

[59]  HR3*-  12  14  16  18  20 

[60]  'RANGE  INTERVAL  FOR  SIMULATION:  < 

[61]  'MINIMUM  RANGE  (IN  METERS):  ' ,9RMIN 

[62]  'MAXIMUM  RANGE  (IN  METERS): 

' , VRMIN+ (NSTEPSxRESO ) + ( HNSTEPS x HRESO ) 

[63]  1  SEA  STATE:' 

[64]  SEA 

[65]  o  'SEA'  IS  ROUGHNESS  PARAMETER 

[66]  'FREQUENCY  (IN  GHZ):' 

[67]  F 

[68]  'INITIAL  RANGE  (IN  KILOMETERS):' 

[69]  fll+1000 

[70]  'TRANSMITTER  HEIGHT  (IN  METERS )- ACTU AL :  ' 

[71]  HAT 

[72]  'RECEIVER  ANTENNAS  HEIGHTS  (IN  METERS):' 

[73]  HR3UHR31 

[74]  n  DISPLAING  THE  INITIAL  SIMULATION  DATA 

♦[75]  *0 
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** ST ART PRO 

[0]  F0  STARTPRO  HR ; UI0 ;HA ;R ; LAMBDA ;ET1 ; ET2 ;ET3 ; ETn ; ETS ; ET1 2 ; ET1 3 
£114 ; ET1 5 ; SI ; SI 12 ; ST  13 ; SI  14 ; SI 15;fll2;£13;£l4;#i 5 ;0K12 
OK13  ;OX14  tOKlS  ;  TEMP;  TEMPI  ;TEMP2  ;SH;RD1  ,  °D2  ;RD3 ;P.DH  ; 

RD 5 ; DELTA  1 ; DELTA2 ; DELTA3 ; DELTAS ; DELTAS •  VET A 1 ; TETA 2 ; 
TETA3 ; TETA'X ; TETA 5; XII; XI 2 ; XI 3 ; XI 4 ; XI 5 ; KAO  1 ; RHO 2 ; RHO 3 : 
RHO'* ;RH05 

Cl]  a  MAIN  CALCULATION  AND  COMPUTATION  FUNCTION 

[2]  010+1 

C 3 3  R+RMIN+(HRESOx\HNSTEPS) ,  (.HRESO*HNSTEPS)+RESO*\NSTEPS 

[4]  a  <R'  IS  S4M7S  VECTOR  ( '  HNSTEPS '  STEPS  OF  'HRESO'  METER 

[5]  a  EACH,  PLUS  ' NSTEPS '  STEPS  OF  <RESO<  METER  EACH), 

C6 ]  a  STARTING  AT  ' RMIN' 

C  7 ]  HA+MINHEIGT+HEICTRESx i HEIGTST 

C8]  a  'HA'  IS  TRANSMITTER  HEIGHT  VECTOR  ('HEIGTST'  STEPS  OF 
[9]  a  ' HEIGTRES 1  METER  EACH ) 

CIO]  GRH+ (HEIGTST, ( NSTEPS + HNSTEPS ) ) p 1 

Cll]  a  'CRH'  IS  MATRIX  WHICH  SHOWS  ALL  THE  LOCATIONS  OF 

C12]  a  CANDIDATE  RANGE  AND  HEIGHT  COMBINATIONS 

C13]  LAMBDA+0 .3*F0 

C  1  4  ]  A 

C15]  SD1+(((M-HSC1]  )*2  ) 0  .  +R*2  )*0 . 5 

C 16 ]  RD2+(((HA-HR(2) )*2 ) « . +£* 2 )*0 . 5 

C17]  RD3+(((HA-HRC3) )*2 )« . +R*2 )*0 . 5 

C18]  i?I54+(((ffil-ffi?C4]  )*2  )°  .  +R*2  )*0 . 5 

C19]  SD5*(((SA-flSC5] )*2 )« . +R*2 )*0 . 5 

C20]  a  'RD'  IS  THE  DIRECT  DISTANCE  BETWEEN  TRANSMITTER 

C21]  a  AND  RECEIVING  ANTENNAS  ( RESPECTIVELY ) 

C22]  DELTAl+(2*HAxHRlll )».*R 

C23]  DELTA2+(2*HA*HRt2] )°.*R 

C24]  DELTA3+(2*HAxHRC3] )».*R 

C 2 5 ]  DELTAm-(2xHAxHR[H]  )o.*R 

C26]  DELTA5*(2xRAxHR:5] )• .*R 

C27]  a  'DELTA'  IS  THE  DIFFERENCE  BETWEEN  THE  DIRECT 

C28]  a  PATH  AND  THE  REFLECTED  PATH  FOR  EACH  RX . 

C29]  o  ANTENNA  ( RESPECTIVELY ) 

C30]  TST>ll*"lo(^(<t>p/?Pl)p \HR(\)-HA  )*RD  1 

C3 1]  r£M2-*-”lo($(4>p/?D2)p  \HR(2)-HA)*RD2 

C32]  TETA3+~la (§($pRD3 )p I HR131-HA )*RD3 

C33]  rE2*A4+'lo(^(4)pSS4)p  |MC4]-Hi4)+SD4 

C34]  r£W5*’lo(^(<))p/?05)p!SflC5]-ffA)tSD5 

C35]  o  'TETA'  IS  THE  OBSERVATION  ANGLE  BETWEEN  THE  RX. 

C36]  a  ANTENNA  AND  THE  TRANSMITTER  (FROM  THE  LATEST 

C37]  a  POINT  OF  VIEW)-  FOR  EACH  RX.  ANT.  (RESPECTIVELY) 

C38]  XI\+(*i($pRDl)pRAl\)+HA)*RDl+DELTAl 

C39]  XI2+(Sl(<bpRD2)pHRl21+HA)*RD2+DELTA2 

C40]  XT3+($($pRD3 )pHRL3)+HA)*RD3+DELTA3 

C 4 1  ]  XI 4* ( H(4>pRDu  )pHR(u)+HA)  *RD'*  +DELTA 4 

C42]  XIS+(ts>(<PpRDS)pHRCS]+HA)*RDS+DELTA5 

C43]  a  <XI<  IS  THE  SINE  OF  THE  GRAZING  ANGLE  FOR  EACH 

C 4 4 ]  a  RECEIVING  ANTENNA  (RESPECTIVELY) 
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♦[45] 
♦  [46  ] 
♦[47] 
♦[48] 
•*■[49] 
♦[50] 
♦[51] 
♦[52] 

[53] 

[54] 
♦[55] 

[56] 

[57] 

[58] 
♦[59] 

[60] 

[61] 

[62] 

♦[63] 

[64] 

[65] 

[66] 
♦[67] 

[68] 

[69] 

[70] 
♦[71] 

[72] 

[73] 

[74] 
♦[75] 

[76] 

[77] 

[78] 
♦[79] 

[80] 

[81] 

[82] 

♦[83] 

[84] 

[85] 

[86] 
[87  J 
[88] 

[89] 

[90] 

[91] 

[92] 

[93] 

[94] 


*(SEA=  0  )/S 0 
♦(SEA=1 )/Sl 
♦(SE4= 2 )/S 2 
♦(S£A=  3  )/S 3 
*(SEA= 4 )/S4 
♦(SE4=5 )/SS 
*(SEA=6)/S6 
■*(SEA  =  7)/S7 


S0:SH*SHL*SHH*0 

+CONT 

51  :Sff<-0 .0325 
SHL+ 0 
SHH* 0.065 
+CONT 

S2:Sff+0.1375 
SSL* 0.065 
SHH* 0.21 
■+CONT 

S3  :  SH*0 . 26  5 
SHL+0 . 21 
SHH+0 . 32 
+CONT 
S4:SH+0 .43 
S55+0.32 
SHH*0 . 54 
+CONT 
S5 :SH*0 . 7 
SSL* 0.54 
Sffff+0.86 
*CONT 
S6:S5^1.08 


SHL*0 . 86 
Sfffl+1 . 3 
-+CONT 
S7 : SH+1 . 7  5 


SHL*1 . 3 


Sff5^2 . 2 
+CONT 

o  'SM'  AND 
o  (.RESPECTIVELY)  RMS  WAVES 
CONT : 


THE  LOWEST  AND  HIGHEST 
HEIGHT  FOR  EACH  SEA  STATE 


RH01+*~2x((o2 )*SH*XI1*LAMBDA )*2 
RH02+*~2 x ( (o2 )*SH*XI2*LAMBDA )*2 
RH03+*~2x ((o2)xSHxXI3*LAMBDA)*2 
RH0Li**~2x  (  (o2  )xSHxXln*LAHBDA  )*2 
RHO 5+*~2x((o2 )xSHxXI5*LAMBDA )*2 
0  -75  2775  REFLECTION  COEFFICIENT  FOR  EACH 

R  SX.  ANTENNA  (RESPECTIVELY)  -  CALCULATED  AFTER 
a  £MEN£  MODEL  (SEE  7 SIM) 
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[95] 

[96] 

[97] 

[98] 

[99] 
[ICO] 
[101] 
[102] 

[103] 

[104] 

[105] 

[106] 

[107] 

[108] 

[109] 

[110] 
[111] 
[112] 

[113] 

[114] 

[115] 

[116] 

[117] 

[118] 

[119] 

[120] 
[121] 
[122] 

[123] 

[124] 

[125] 

[126] 

[127] 

[128] 

[129] 

[130] 

[131] 

[132] 

[133] 

[134] 

[135] 

[136] 

[137] 
♦[138] 

[139] 

[140] 
♦[141] 
♦[142] 

[143] 

[144] 

[145] 

[146] 
♦[147] 


ETl+RHOl*  (lo(((o2  )x(  (HA*2  )+SB[  l]x«i  ))«.■»  (LAMBDAxR))  )*2 
ET2+F.H02*  (lo(((o2  )*((HA*2  )+SB[2]  *HA  ) )  »  .  i  (LAMBDAxR  )  )  )*2 
ET3+RH03 x(io(((o2)x((B«*2  )+fffi[3]*M ) )° . i (LAMBDAxR )  ))*2 
£T4-7?fl04x(io(  (  (o2)x((M*2)+fifi[i4]*M  )  )<> .  *  (LAMBDAxR  )  )  )*2 
ETS+-RH05x  (lo(((o2  )x  (  (M»2)tffl[5]*M  ) )°  .  *  (LAMBDAxR  ) )  )*2 
o  ' ET 1 1  +  1  £T5 1  ABB  MULTIPATH  ATTENUATION  FACTORS  FOR 

a  EACH  RECEIVING  ANTENNA,  RESPECTIVELY  ( MATRIX ) 
ET12*-ET1*ET2+1E~13 
ET13+ET\*ET3+1E~13 
£T14+Eri+E2’4+1£-13 
£Tl5+£2’lt£Z’5  +  lB-13 

a  'Bn 2 ' * 'ET1S '  ABB  RELATIONS  BETWEEN  THE  MULTIPATH 
a  FACTORS  FOR  EACH  PAIR  OF  RX .  ANTENNAS,  RESPECTIVELY 
A  (MATRIX) 

RCV-.SI+SIM 

a  'SI'  ARE  RECEIVED  SICNAL  LEVELS  FOR  ALL  THE  RECEIVING 
a  ANTENNAS  ( VECTOR ) 

REL: 

SJ12+(pCBB)pSJ[l]+SJ[2]+lB“l3 
SJ 1 3-*-  (  pCBB  )pSJ[l]  +5/  [3  ]  +lB”l  3 
SI14^(pCBB)pSI[l]*S/[4]+lB"l3 
SJl5*(pCBB]pSI[l]+SI[S]+lB"l3 
a  ' SJ1 2 1  * ' SI IS '  ARE  RELATIONS  BETWEEN  RECEIVED  SIGNALS 
a  PAIRS  (MATRICES ) 

COMP : 

R12+SI12*ET12+1E~13 
R13+SI13*ET13+1E~13 
B14+SJ14+£2’14  +  1B“13 
Bis^sns+Bns+iB'is 

a  'B12'*'B15'  ARE  RATIO  MATRICES  BETWEEN  ACTUAL 

a  AND  THEORETICAL  VALUES 

OK12*-(R122C  )A  (R12<,D  ) 

OK13*(B132C)A(B13SB) 

0B14+(B14SC)a(R14sI3  ) 

OB15+(B152:C)a(B15SD) 

a  'OK12' *'OK15'  ARE  LOCATIONS  OF  RECEIVED  SIGNALS 
a  WHICH  ARE  EQUAL  TO  THE  PRE-CALCULATED  VALUES 
a  (ACTUAL  RELATIONS  EQUAL  TO  THE  THEORETICAL  ONES) 
a  (MATRICES) 

MATCOMP: 

CBB+OK12aOB13aOB14a0/C15 

a  'CRH'  ARE  LOCATIONS  WHICH  ARE  EQUAL  ALSO  BETWEEN 
a  ANTENNAS  PAIRS 
♦(1=+/ ,CRH) /RHFIN 

a  CHECK  IF  THERE  IS  A  SINGLE  COMBINATION  OF  RANGE 
a  AND  HEIGHT 
■*(FLAC-\  ) /RHFIN 
♦BCV 

a  IF  MORE  THAN  ONE  COMBINATION ,  CO  TO  ANOTHER  RECEIVING 

a  SESSION  (AS  LONG  AS  THE  RANGE  IS  RICHER  THAN  THE 

a  MINIMUM  RANGE).  IN  CASE  OF  A  SINGLE  COMBINATION  CO  TO 

a  TRANSLATE  INTO  REAL  READING  (RANGE  AND  HEIGHT) 

RHFIN -.*0 
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**SIM 

[Cl  A*SIM;LAMBDA:A  1  ;A2  ;A3  ;  44  ;  A5  ;  SXI ;  SHR ;  SRHO 

[1]  o  RECEIVED  SIGNAL  AMPLITUDES  SIMULATION  ( AT  EACH 

[2]  e  RECEIVING  ANTENNA ) ,  INCLUDING  TRANSMITTER 

[3]  o  MOVEMENT  INFLUENCE 

[u]  LAMBDA*0.3*F 

[5]  FLAC*Q 

[6]  SXI*-~  3  o  (  (  |  HR3  -HAT )  *RI ) 

[7]  o  'SXI'  IS  THE  CRAZING  ANCLE 

-*■[8]  +(SEA-0)/LABl 

[9]  SHR*SHL+(1?  (SHH-SHL  )xlOOO)'MOOO 

CIO]  o  'SHR'  IS  THE  RMS  WAVES  HEIGHT  FOR  SIMULATION 

->-[11]  *LAB  2 

C 12]  LAB1:SHR*0 

[13]  LAB2:SRHO**~2*( ( (o2 )*SHR*SXI* LAMBDA  )  )*2 

[14]  o 

Cl 5]  a  THE  SIMULATION  USES  AMENT  MODEL  FOR  CALCULATING 

[16]  p  THE  REFLECTION  COEFFICIENT  (p): 

[17]  p 

[18]  p 

[19]  p  p  =EXP.  (-2*(2*PIxH  *SIN . (XI)*LAMDA)*2) 

[20]  p  S  RMS 

[21]  p 

[22]  p 

[23] 

Al+SRHO  [  1]  x  (FFPxUx  (lo(((o2  )x  (HR3  [1]  +HAT  )*HATU  (LAMEDAxRI )  )  )*  2  )*  ( 
04  )y.RI* 2 

[24] 

A2+SRH0121  *  (FFPxux  (lo(((o2)x  (HR3  [2]  +HAT)*HAT )  *  (LAMBDA*RI )  )  )*2  )*  ( 

04 )xflj* 2 

[25] 

A3*SRH0  [  3  ]  x  (FFPxu  x  (io(((o2  )x  (HR3  [3]  +HAT)*HAT)±(.LAMBDA*RI ) )  )*2  )*( 
04  )x/?J*2 

C  26  ] 

A4«-Sfltf0[4]x(£7?px4x(io(  ((o2  )*(HR3('i1+HAT)xHAT)*(LAMBDAxRI)  )  )*2  )*  ( 
04 )xRI*2 

[27] 

A  5+SRHO  [  5]  x  (FFPxux  (lo(((o2  )*  (HR3  (51 +HAT  )*HAT  )*  (LAMBDA*RI )  )*2)*( 
o4 )xRI*2 

C 2 8 ]  P  'SPflO'  IS  THE  'REFLECTION  COEFFICIENT'  (p)  FOP 

[29]  p  EACH  RECEIVING  ANTENNA  ( RESPECTIVELY ) 

C  3  0]  A<-i41,A2fA3,A4,A5 

[31]  p  'A' IS  A  VECTOR  OF  RECEIVED  AMPLITUDES  (.IN  EACH 

[32]  R  ANTENNA -RESPECTIVELY ) 

[33]  RI+RI-SIMSTEP 

[34]  p  UPDATING  THE  TRANSMITTER  RANGE  (IN  'SIMSTEP'  METER  STEPS ) 

-[35]  - (RI2RMINS ) /RET 

[36]  '  ' 

[37]  'END  OF  SIMULATION  SESSION’.  < 

[38]  FLAC-l 

[3  9]  p  RANCE  HAD  NOT  BEEN  FOUND  YET,  BUT  IS  LFSS  THAN  THE 

[40]  p  MINIMUM  RANGE,  OR  END  OP  SIMULATION  SESSION 

-[41]  RET:* 0 
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«*RHFIND 

Co:  RHFIND ;  INDH\  INDR;MATHR;  INDVEC;  LOC:  VECCRH 

Cl]  o  FINDING  AND  DISPLAYING  THE  RANGE  AND  HEIGHT 

[2]  a  OF  THE  TRANSMITTER 

C  3 ]  INDH* 100000 XMINHEIGT+HEICTRES * i HEIGTST 

Cu]  a  1 INDH '  IS  AN  HEIGHT  INDEX 

C  5 ]  INDR+RMIN + ( HRESO x i HN STEPS ) , ( HRESO x HN STEPS )+(RESO*\ N STEPS ) 

C6]  a  'INDR'  IS  A  RANGE  INDEX 

[  7  ]  MATHR+INDH  •> .  +INDR 

[8]  a  1  MAT HR 1  IS  A  RANCE  /  HEIGHT  INDEX  MATRIX 

[9]  INDVEC*, MATBR 

[10]  FIND : 

[11]  V EC CRH*. CRH 

[12]  a  CHANGING  MATRICES  INTO  VECTORS 

[13]  LOC*VECGRH / INDVEC 

[14]  a  SELECTING  THE  CORRECT  LOCATION  OF  THE  TRANSMITTER 

[15]  a  ('LOC'  IS  THE  TX.  LOCATION ) 

[16]  H <-100000  I  LOC 

[17]  ff<-LLOC*100000 

[18]  a  TRANSLATING  INTO  REAL  RANGE  ('R')  AND  HEIGHT  Off') 

[19]  " 

[20]  '  ' 

[21]  'TRANSMITTER  RANCE  (IN  KILOMETERS):' 

[22]  fl+1000 

[23]  'TRANSMITTER  HEIGHT  (IN  METERS):' 

[24]  H 

[25]  '• 

[26]  '  ' 

[27]  'UPDATED  ACTUAL  RANGE  (IN  KILOMETERS):' 

C261  (SIMSTEP*10Q0 )+RI*1000 
♦[29]  ♦0 
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